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I. INTRODUCTION 


The Great Basin Restoration Initiative (GBRD is a multi-state, multi-agency, state, local, 
and federal initiative focused on “restoration” of the sage/pinion/juniper biome within 75 
million acres of the Great Basin. To advance restoration efforts, ecological planning 
boundaries were created and multiple-use management goals, including those pertaining 
to cultural resource management and protection were, or are in the process of, being 
developed. Many of the management goals of the GBRI can be accomplished if cultural 
resources (significant historic, prehistoric, and ethnohistoric sites and localities) can be 
managed in a more efficient manner. 


One way to facilitate management and planning is to develop and test cultural resource 
distribution models that predict site density and distribution for planning purposes. Over 
the last several years, predictive models have been generated for relatively large 
hydrologic basins (Railroad Valley and Pine Valley) in Nevada (Drews et al. 2002; 
Zeanah 1998) and Utah (Zeanah 2001). These models are based on relatively fine-grained 
analyses of landscape, soils and geomorphology and they predict cultural trends that 
appear to be valid within their respective hydrologic basin. 


While landforms and vegetation classes are relatively consistent across the Great Basin, 

orthographic effects of bounding mountain ranges create microclimates within each 

hydrographic basin so that vegetation and landform mosaics are not always comparable 

across broad areas of the landscape. The challenge for the GBRI cultural resources model 

is to test a larger area, coarser, landscape level modeling. Based on the model, areas & 
within the GBRI area can be more effectively managed to ensure efficient use of a BLM 

district’s resources while furthering the goals of the GBRI. The model should be a basis 

for understanding history and prehistory of the GBRI landscape and how humans have 

positioned themselves on the landscape over time. 


With that task in mind, an extensive project boundary consisting of 12 major 
hydrographic basins covering 20,533,700 acres within Nevada, Utah, and Idaho was 
chosen for study. The area includes environments typical of Nevada basins, the Snake 
River plateau, and the Great Salt Lake basin. Model results are useful to managers at the 
Bureau of Land Management, Elko and Ely field offices in Nevada; the Salt Lake, 
Filmore and Cedar City Field Offices in Utah; and the Owyhee and Jarbidge field offices 
in Idaho. 


A Perspective on Study Goals 


The goals of this study are, to be frank, more managerial than anthropological. Our 
charge was to examine where archaeology is most likely to be found, and when found, 
where it was most likely to be an impediment to Great Basin Restoration Initiative land 
use goals. The management orientation of this study does not mean there is no 
component of science to it. Our ability to formulate reasonable hypotheses about historic 
and prehistoric settlement patterns — to devise simple initial models at all — derives from 
research by ourselves but especially by others. We make no claim to be doing “deep 
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science” in the hypothetico-deductive mode in this study. Nevertheless, we have striven 
to make the results and information useful to those interested in more focused work of 
that sort. 


The sheer size of the study area itself precludes all but the most cursory of 
scientific,deductive, model-building on human behavior. Typically, one requires very 
fine-grained information for particular time periods to create an effective deductive study. 
For instance, Zeanah (1995) used an area approximately 5% the size of this study area for 
his analysis of prehistoric foraging patterns and the resulting archaeological record in 
Churchill County, Nevada. As we discuss below, the study has utility for management 
and for researchers, but in different ways. 


Management and Implementation Goals _ 


A central concem with most models is how managers will interpret and implement the 
results. The goal of this project is xot to create a lock step management document (€.g., 
prescribed treatment within specific areas), but rather as a planning tool for BLM 
managers, biologists, and cultural resource specialists to evaluate potential conflicts as 
they work within the GBRI. From this general perspective, several research questions 
can be generated. 


e Which landscape factors are the best predictors of cultural resource location? 

e What management characteristics in terms of National Register status or Cultural 
Resource Use Allocations do resources have, and how are they distributed? (cf. 
BLM Manual Section 8110, “Identifying Cultural Resources”) 

e Do sensitivity boundaries relate to criteria that are readily observable in the field 
and can they be identified through simple overlay of available data? Are more 
complex analyses required for the model to be effective? 

e Is the planning model a useful tool to aid in the identification of areas where 
imminent threats from natural or human caused actions may cause deterioration of 
significant cultural resources? | 

e Defining limits of knowledge and areas of further information needs. 

e Procedural recommendations concerning subsequent data gathering, testing and 
strengthening of the model. 


The work also provides a chance to contrast the management outcomes of broad-scale, 
inductive models with detailed deductive models based on optimal foraging theory. 
Railroad Valley, Pine Valley, and the Dugway Proving Grounds area all lay within or 
nearby the study area. These three areas are or have been examined with detailed forager 
behavior models. The anthropological models have resulted in a consideration of 
management plans and needs within respective study areas. Those models provide an 
informative contrast to this large-scale model. 
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Anthropological Goals 


The large, and diverse, scientific literature on human settlement patterns is the 
anthropological context of this study. Although the study goals are, in essence, the 
discovery of correlation and not its explanation, settlement pattern is necessarily an 
outcome of the study too. This study undoubtedly raises more questions than it attempts 
to answer. That is one of the shortcomings of correlation approaches in general. 
Nevertheless, our intent is that the questions that it raises are themselves useful scientific 
leads for further research. 


Some of the questions that a study such as this leaves unanswered but tantalizingly 
available for speculation include: Are there discernable relationships between site 
assemblages and the landscape? Springs and other perennial or predictable seasonal water 
sources are thought to be attractants of prehistoric use. Is this really the case? What about 
the potential for buried sites; is this usually greater near springs? Likewise, is there a 
discernable pattern of early sites along Late Pleistocene/Early Holocene lakeshore 
margins? Are potential wetland environments good predictors of sites density? 


Beyond settlement pattern studies themselves, there are questions of change through time 
in the prehistoric archaeological record that the study examines in broad view. 
Subsistence change was one such question that we hoped to address. The occurrence of 
pinyon pine in the study area and its role as a dietary staple generates several issues of 
interest. Over time, pinyon has general spread from south to north, fluctuating up and 
down in elevation in response to temperature and precipitation regimes (Grayson 1983). 
Pinyon may have never been present in the northernmost portion of the study area. 


We assume that there is a detectable archaeological signature for pinyon exploitation. In 
many parts of the basin, that signature is rock rings and groundstone implements, within 
proximity to the current pinyon-juniper zone (see Thomas and Bettinger 197 6:272). Sites 
containing those features may indicate the overall range of pinyon through time, in terms 
of its expansion from the south, as well as localized elevational expansion and 
contraction. 


This assumption may be faulty since groundstone could have been used to process any 
number of seed resources and all vegetation expands and contracts with climatic 
variation. Site density alone may be a better measure of exploitable resource zones. 


Several antelope traps are known from this part of the Great Basin. The west central 
portion of the study area is characterized by high mountains that collect moisture during 
the winter and with relatively low, open valleys that come into production during early 
spring and are moistened by runoff late into the summer. Valleys provide ample forbs 
and browse. Open juniper woodlands in the foothills provide access to construction 
materials for drive fences. Analysis of landscape in the vicinity of known antelope drives 
may serve to develop a testable hypothesis for site location. 








Potential grasslands in the Snake River uplands may have provided prime bison habitat. 
A correlation between grasslands and Northern Side-notched projectile points might 
suggest big game hunting. 


The distribution of Fremont sites across the landscape may reveal land use patterning and 
contingencies for site location relative to productive agricultural lands or specific 
resource procurement areas. Likewise the distribution of Late Prehistoric and pre-contact 
projectile points may identify Numic progressions from south to north across the project 
area. 


Because of the size of the project area and the goals of the GBRI, the cultural resources 
models (one for prehistoric, another for historic period) is not spatially or temporally 
fine-grained. We used hundreds of spatial units for soils, vegetation, and topography. 
Yet, the size of the study area was so large that even the smallest spatial units are equal in 
size to the largest spatial units in for example, the Pine Valley study. Large area, 
landscape-level datasets were tested as predictors of cultural resource distribution and 
significance. In a sense, the study tests both the correlation itself and the methodology: 
the validity of developing models over such a Jarge area using using spatially and 
categorically coarse datasets. 


Report Overview 


The study report has a simple structure. We first describe the project setting, including its 
natural and historical contexts (Chapter II). Next, we discuss how analytical units were 
divided out of the study area as a whole (Chapter III). Chapter I'V is an exposition of the 
study methods. Chapter V presents the bulk of the study results, followed by discussion 
and closing comments (Chapter VJ). 


Il. MODEL BUILDING: A BACKGROUND e 
x 


Although the study goals are not anthropological entirely, our ideas about how best to 
produce a “forecast” of where one was most likely to find archaeological sites derived 
from the anthropology of the Great Basin. Anthropology — the study of man — has always 
been closely tied to archaeology in the Great Basin. Our approach to model-building 
drew from the seminal work of Julian Steward, as did the work of many, many, other 
archaeologists in the region. 


In this chapter, we present brief discussions of these major models. Relatively little of 
what we are doing here is new, novel, or untried. So, these short summaries inform as to 
why we took particular decisions in this research. 


Environmental Models 


Julian Steward laid the groundwork for much of the research done in the Great Basin 

subsequent to the publication of his seminal work Basin-Plateau Aboriginal 

Sociopolitical Groups (1938). Researchers used the results of Steward’s ethnographic 

reports as springboards to study the environment as a limiting factor in the level of 

cultural complexity attained by people in semi-arid landscapes. Limited resources in 

these areas would force the inhabitants to spend the majority of their time and effort 

procuring food and producing the technology required to aid in these tasks. Subsistence 

and settlement patterns could then be explained and explored in these terms. Steward’s 

information provided many research topics and continues to be a valuable source of @ 
information for archaeologists. 


Archaeologists quickly picked up from Steward the importance of pinyon nuts as a staple 
resource. Several studies examined whether the ethnographic reliance on pinyon had 
antiquity (Thomas 1971, 1973; Thomas and Bettinger 197 6). Central Nevada studies 
showed that pinyon was, indeed, one of the long-term staples. This led to a focus on 
pinyon as the determinant of prehistoric settlement pattern. 


More generally, Steward’s work created a family of models that Wilde (1994) describes 
as paleoecological models, including several relating to the Great Basin. These include: 


@ Steward’s ethnographic model discussed above which argued that the Great Basin 
had a “socially fragmenting effect upon its prehistoric inhabitants”; 


e Jenning’s “Desert Culture” model based explicitly on Steward’s work, originally 
set out to account for the record at Danger Cave in northwestern Utah and 
proposed a cultural ecological model in which a stable settlement and subsistence 
pattern was evidenced for the past 10,000 years; 


© The Wamer Valley model as delineated by Weide (1968) which is a lake and 
marsh-oriented pattern, but with increased reliance on upland faunal resources, 
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® O’Connell’s Surprise Valley model with early and later variants. The early period 
(6500-4500 B.P.) was based primarily on marsh and grassland resources 
supplemented by upland animals during the winter; and, 


© The Steens Mountain Model which shows an inverse relationship between site 
frequency and site size, which suggests high resource productivity, allowed 
population aggregations (Wilde 1994:97-102). 


For a period of about 15 to 20 years, the pinyon-centric model of aboriginal settlement 
was truly the dominant paleoecological model. Because pine nuts are easily harvested 
and give a rich return, their absence can be predicted up to two years in advance, and 
their presence at least predictable in location, they were given primacy in many of the 
models of Great Basin prehistory. However, further archaeological work showed that the 
pinyon-centric model of aboriginal subsistence and settlement was too narrow. Other 
natural settings in the Great Basin, especially wetland and lacustrine environments, have 
long and rich archaeological records too. Resources in these settings are not so easily 
understood as pinyon nuts. For instance, why would one gather cattail pollen instead of 
harvesting pinyon nuts? Was this an alternative subsistence strategy or equal to “king 
pinyon”? Exploration of these questions brought anthropologists and archaeologists in to 
a consideration of caloric maximization in patchy environments: optimal foraging theory. 
Optimizing theory attempts to understand, and thus predict, the choices that a rational 
forager will make. 


Optimizing models have been very successful as a deductive form of environmental 
model. Many studies in the relatively stark Great Basin have used models of what 
foragers should have done as rational behavior. These studies then examine whether the 
archaeological record matches the predicted behavior. Generally, such studies have been 
successful over areas of about half a million acres, such as a typical basin and range 
valley (Bonstead 2000; Connolly 1999; Gehr 1980; Jones et al. 2002; Mehringer 1986; 
Nials 1999, 2000; Pendleton 1979; Pettigrew 1984; Pinson 1999; Thomas 1971). 


In the central Great Basin, the primary GIS optimal foraging models for the Great Basin 
have been proposed by Zeanah et al. (1995), Zeanah (in press), and Raven and Eiston 
(1989). Beck and Jones (2000) provide a thoughtful overview of how these efforts fit 
within contemporary regional research directions. Optimal foraging approaches are not 
without problems. One of the main criticisms of the use of these models is that they are 
not easily replicated. Though they go far in description, they offer little in explanation 
outside of resource return rates in the form of calories expended and/or gathered per hour. 
Optimal models provide detailed formulas for energy return rates, but do not account for 
resources used in other contexts such as medicine, ritual, fuel, or shelter. 


Overall, then, the history of inquiry in Great Basin archaeology has gone from informal 
paleoenvironmental models to ever more detailed and quantitative approaches. The latter 
methodologies, especially optimal foraging, provide numeric baselines from which to 
understand prehistoric settlement patterns. ; 
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Geomorphic Approaches 


A geomorphic site preservation approach has been applied to archaeological studies in 
the Great Basin. This model developed and refined by Nials (1999; 2000) uses 
geomorphic principles to identify areas likely to retain in situ cultural materials. Suitable 
locations include those lying on and adjacent to: 


« Late shorelines of pluvial lakes, including dunes contemporary with late pluvial 
lake shorelines; 


e Distributary drainages entering open basins; 


© Upland valley bottoms where stream gradient locally flattens out and the valley 
widens; 


e Near springs active at the appropriate times and; 
e Rockshelter and caves (Nials 1999). 


This approach contains some tautological assumptions pointed out by Nials (personal 
communication; 2002) that make it problematic for elucidating patterns in the 
archaeological record. The geomorphic model promotes the survey of landforms that are 
favorable for, and have a high probability of, containing intact sites. In other words, well- 
preserved sites are looked for in the exact environments in which they should be found. 
Whitley (2000) notes that cause and effect in the record become difficult to discern: 


For instance, correlating 97% of sites with floodplains is meaningless if 
97% of the survey areas from which the data is derived occur on 
floodplains. . . it is assumed that geomorphological setting was a constraint 
on site locations for instance, yet it is rarely clear how important certain 
geologic structures were in comparison with the relationship to a permanent 
source of water. Secondly, it is unclear whether it is the geomorphological 
setting or the distance to water is important, if there is already a spatial 
correlation between the two. . .Whitley (2000:27). 


Approaches to Model Formation 


The analysis done for the GBRI project does not attempt to falsify other models but 
points to the fact that they may not be the best approaches for GBRI. Selectionist models 
and others that attempt to explain human behavior in terms of natural selection are based 
on biological principles of animal behavior. In many cases, applying that theoretical 
approach to model human behavior does not provide adequate explanations for social 
components of the system, though they do provide general descriptive frameworks useful 
in explaining optimal utilization within broad environments. Perfect information about 
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people’s environment is rarely available, “which means that they never really forage 
optimally, but base decisions on their best guesses” (Kelly 1995: 100). Optimal models 
work well in small, delimited areas, but would not be practical to derive for the 20 
million acres under consideration for the GBRI project area, which covers parts of three 
states in different topographic and environmental settings. 


Because the units of analyses are hydrographic units (HUC) explored in a Geographic 
Information System (GIS) environment, the next section provides a very basic discussion 
of the premises of GIS. GIS will be used in the plural when referred to in a general sense, 
as there are many GIS programs, and in the singular when used in reference to the results 
of this particular project, as the final models were built using one particular product, 
ArcView® version 3.3. 
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Ii. STUDY AREA DEFINITION 


The study area occupies the north-central portion of the Great Basin physiographic 
province and straddles the states of Nevada, Utah, and Idaho (Figure 3.1). After accessing 
data quality and completeness of GIS and database files for the region, a final 
configuration, consisting of 12 hydrographic basins was decided upon (Table 3.1) 
Hydrographic unit boundaries are derived from United States Geological Survey data, 
compiled from 1:250,000 base maps in 1973. The study area encompasses an excellent 
cross-section of the Great Basin biome, though a small portion of the area lies outside of 
the Great Basin as a hydrographic entity (Figure 3.2). For the sake of simplicity, we have 
used the common abbreviation of “HUC” (hydrologic unit catalog) in much of the 
following reporting to mean an individual hydrologic unit catalog item, 1.¢., a 
hydrographic basin. 


Our initial attempt at model development proposed that we employ broad, landscape 
level criteria to model the probability of encountering cultural resources over the broad 
reach of the entire study area. While grossly similar in terms of very general climate, 
vegetation, and topography, considerable variation is apparent as one moves from west to 
east and south to north across the study area. When contrasted by size, landform, and 
hydrologic regimes, the valleys in eastern Nevada bear little resemblance to either the 
Great Salt Lake basin, or the Upper Snake area. In order to maintain the landscape level 
approach but constrain environmental factors, the sub-basins within the study area were 
chosen as the analytical units rather than the larger study area as a whole. 


Variations in topography define three hydrographic sub-regions within the study area; the 
Central Nevada Desert, Great Salt Lake and the Upper Snake. Each sub-region consists 
of a number of smaller hydrographic units, with boundaries based upon drainage patterns. 
Gross similarities in topography and drainage pattern and vegetation allowed for analysis 
of the Upper Snake and Great Salt Lake sub-regions as single units. 


Internally drained basins and surrounding mountain ranges comprise the Central Nevada 
Desert sub-region. While grossly similar, each hydrographical basin within this sub- 
region contains significant variation in physiography that could effect human adaptation. 
Two hydrographic units in Nevada, Long/Ruby Valley and Spring/Steptoe Valley were 
treated separately so that we might better highlight variation between each area, and 
directly compare the results of this landscape modeling effort with more fine-grained 
modeling recently conducted in the Nevada basins of Pine Valley and Railroad Valley. 
The Pilot-Thousand Springs hydrographic unit within the Great Salt Lake sub-region 
showed a closer affinity to Central Nevada units and was likewise analyzed as a distinct 
entity. 
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Table 3.1 
Great Basin Restoration Initiative Analytic Units 
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Figure 3.2 Great Basin Extent 
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IV. METHODS AND PROCEDURES 


Predictive cultural resource models are ‘‘a simplified set of testable hypotheses, based either 
on behavioral assumptions or on empirical correlations, which at a minimum attempts to 
predict the loci of past human activities resulting in the deposition of artifacts or alteration of 
the landscape” (Kohler 1988:33). Based upon their accumulated experience, most 
archaeologists could, on cursory review of a topographic map, accurately predict with 50% to 
80% accuracy where archaeological sites would most likely occur. Predictive capacity alone, 
however, fails to meet the explanatory capacity of rigorous scientific inquiry. Sites that fall 
outside of the predictive pattern are often of greater interest to archaeologists. To better 
understand and evaluate outliers, one must first have a quantitative means to evaluate those 
sites that fall within a “normal” distribution (Heidelberg (2001:6). 


As many as four approaches have been employed as a means identify patterns within 
predictive layers: inductive, deductive, intersecting, and weighted. The inductive approach 
establishes conclusions based upon data that has already been collected. It is widely used 
because it draws from an accumulation of survey and data collection compiled by various 
agencies and researches within a specific region. Biases are inherent due to variable survey 
strategies, sampling criteria and vagaries in data collection methods. Nonetheless, benefits 
are derived reduced costs from utilizing existing data. 


Deductive pattern identification is derived from data specifically collected for the purpose of 
the study. Sampling strategies are controlled and data collection is consistent throughout the 
model area allowing negative findings to be more readily assessed. Additional background 
layers consisting of regionally specific data on vegetation, elevation, slope, aspect, soils, 
hydrology, and climate can be used to test deductive hypotheses regarding human land use 
decisions. 


An intersecting approach combines deductive or inductive data sets with background layers 
to define probability within each environmental layer. When several probability zones 
overlap, their intersection defines an area of high sensitivity, with fewer overlaps defining 
medium and lower sensitivity zones. 


A significant problem with the intersecting approach is that all variables are considered 
equally. To counter that shortcoming, environmental variables can be weighied so that a 
theme, such as aspect, is considered a lower relative value than say, distance to water. A 
scalar variable may also distinguish relative values within each environmental class. 
Combining intersecting and weighting methods creates an even more robust approach. 


Model Processes 


Determining which environmental and cultural variables and how those variables would be 
analyzed was a major consideration for the development of the planning model. Initial test 
runs with a limited data set from northeastern Nevada utilized chi-square analysis to 
determine the distributional relationship of sites to distinct environmental zones. The process 
required extensive manipulation of tabular and grid data sets then subsequent overlay of 
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predictive themes to produce a generalized sensitivity map. Updates and model testing 
utilizing this method would require continued technical expertise, thus reducing the overall 
utility of the model as a planning tool. A more economic approach to modeling was sought, 
in which new data could be easily input and new models generated in response to additional 
information. 


A weights-of-evidence software package, Spatial Data Modeler (Kemp et al. 1999), was 
recently developed to run with the ArcView® Spatial Analyst extension (ESRI, Redlands, 
CA). It integrates a number of not-dissimilar steps used in the initial chi-square analysis and 
showed promise as a user-friendly, programmatic approach to developing a predictive model. 
To test its reliability and to better understand the modeling program, we contrasted the 
weights-of evidence results with a cell based chi-square analysis. 


Spatial Data Modeler (SDM) is an Arc View® extension developed by the United States 
Geological Survey for mineral exploration purposes. SDM has several options for exploring 
data including: weights of evidence (WofE), logistical regression, fuzzy logic, and neural 
networks. Weights of evidence is particularly useful in predicting mineral deposits based 
upon the location of known resources and archaeologists have successfully applied the 
application to predict the probability of site locations. 


Weights-of-evidence is a discrete multivariate method originally developed in a nonspatial 
context for combining a number of medical symptoms to predict disease (Bonham-Carter 
1994; Xu et al. 1992). “In this situation, the response variable (presence/absence of disease) 
is binary and the predictor variables are also of the presence/absence type.” (Bonham-Carter 
1994:1). Assuming that the variables are not dependent, data sets are combined to give the 
posterior probability to each cell for each unique binary combination. Bonham-Carter (1998) 
explains this idea with the following example: 


If one wished to predict the likelihood of rain for a given day in an area that 
receives an average of 80 days of rain a year, a sound estimate of the prior 
probability of rain would be the ratio 80/365, This initial measure of probability 
can then be modified using other pieces of information to determine the 
probability that it will rain in a particular month depending on the month, the 
focation of the jet stream, or any other factors. The factors determining the 
probability of ram will vary with the time of year and can be figured into the 
equation to produce a model that will answer: “what is the probability that it will 
rain tomorrow?” (Bonham-Carter 1998:302-303). 


Weights-of evidence methods were adapted for use in mineral exploration by overlaying 
geologic and geochemical data sets to predict locations of ore bodies (Bonham-Carter et al. 
1988: Raines 1999), and as a means to predict the location of fossil pack rat middens 
(Mensing et al. 2000). Archaeologist apply this same method in a spatial sense by using 
archaeological sites in an area as training points to create a probability map which aids in the 
prediction of locations likely to contain sites in the area under study. Results can be used for 
numerous purposes but, most recently, have been used by Federal agencies to better manage 
public lands. 





The Bayesian weights-of-evidence approach requires a set of training points, in this case; 
archaeological sites, a set of evidential themes or variables that are assumed to be predictive 
of training point location, and a spatially defined study area. Training points are then 
compared with the evidential themes to calculate a weight assessing the spatial association 
between the points and each class within the theme. A positive weight indicates the class 1s 
present; a negative weight if the class is absent. The strength of a correlation is measured by 
its contrast (W*-W’). Positive contrast values suggest that more training points occur within 
that class than would be expected by chance. Negative contrasts indicate that fewer training 
points within that class than would be expected by chance. The contrast is divided by the 
standard deviation of the contrast values to provide a normalized (Student) contrast for each 
class. 


Positive contrast values are grouped to assess the relative strength of the predictive pattern 
for each class (Table 4.1). Depending upon contrast values, the user determines which 
classes are “inside” (predictive) or “outside” (not predictive) within each evidential theme. 
By determining high or low cutoff points, the user’s decisions directly influence the model 
outcome. In addition, expert opinion can be used to weight an individual class of data 
thought to be intrinsically more important, or to discard contrasts that are artificially high as a 
result of disproportionate unit area to training point values. 


Prior to running the model, the program calculates a prior probability assuming a random 
distribution of sites: 


Prior Probability = Number of training points 
Total of study area units 


Since the training points make up a very small sample of the entire study area, prior 
probability will likely be a number much smaller than the actual density of sites within the 
study area. After weights have been calculated and re-classified into a binary evidentiary 
theme, they are combined to create a response theme that calculates a posterior probability 
for all cells within each unique group of binary combinations. Posterior probabilities that are 
higher than the prior probability suggests a non-random distribution within that intersection 
of evidential themes. 


Information Collection and Evidential Datasets 


Background data used to analyze cultural and landscape features for the planning model were 
acquired from a number of different sources. The challenge with both the cultural and 
landscape data sets was to locate evidential themes that could be applied or adapted to the 
larger study area. In some cases (e.g. geology), consistent data was available for one state, 
but missing from others. Scale was also considered, especially for layers like vegetation, 
where detailed regional coverages lacked comparability between analytic units. 








Table 4.1 
Relative Strength of Contrast Value in Weights of Evidence Analysis 


{after Bonham-Carter; 1994) 


Contrast Value 
P8005 Mildly Predictive 








Cultural Resources and Inventories 


Cultural resource layers compiled for the analysis were derived from a number of different 
sources and required varying degrees of manipulation in order to maximize their utility. 
Idaho and Utah have developed and maintained a geographic information system for cultural 
resources. Both states graciously supplied that information for the project area. Nevada is in 
the process of completing a similar conversion to an electronic archive. As different cultural 
data sets were received, data was merged into a consistent format. All GIS data sets were 
converted from their default projections to a uniform UTM Zone 11, NAD 1927 projection. 


Utah 


Cultural resource shapefiles and resource inventory shapefiles were provided by the Utah 
State Historic Preservation Office. Depending upon relative size of the feature, site and 
inventory locations are displayed as point, line or polygon shapes. For analytical purposes, 
points and lines were buffered to create synthetic polygons and then merged with the 
appropriate (site or inventory) polygon layers to create single polygonal site or inventory 
layers. Attributes for the Utah synthetic shapes included buffered width, area, site or 
inventory number, confidence in plot location, and data entry specifics. Using ArcView 
utilities, a center point was created for each site so that each entity could also be displayed as 
a single point. ; 


Point @ Buffer to 56mradius 9 ——> (e) Buffered point 


——P 
Line “"" Buffer to 15m radius —) Buffered line 


(@) ot —) ~{ ap) as Synthetic Polygon 
Buffered point § Buffered line Polygon 


The Utah site database consisted of a Microsoft Access® database containing Intermountain 
Antiquities Computer System (IMACS) encoded fields. Site numbers in the IMACS database 
allowed the data to be linked to the GIS site shapefiles. 


Idaho 


The Idaho State Historic Preservation Office provided a Microsoft Access® database 
containing UTM coordinates for each site within the project area. Fields pertaining to a range 
of feature types are present in the table structure, and descriptive artifact attributes are 
annotated for each site. A separate table containing SHPO National Register status was 
provided with the site data. Inventory databases with locational information have not been 
compiled for Idaho. 
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Using the Idaho site UTM coordinates, a point theme was created for each site for use in the 
GIS. Attribute tables for the site points contained all tabular data presented in the Idaho 
database. As quarter section data in the inventory database was inconsistent, an attempt to 
determine inventory extent based upon legal descriptions proved futile. Composite legal 
descriptions often produced areas significantly larger than the reported inventory extent, 
making the data unreliable. 


Nevada 


Nevada SHPO maintains site and inventory archives at the Nevada State Museum for its 
northern counties, and at the University of Nevada, Las Vegas, Harry Reed Center for 
Environmental Studies for southern counties. Archival data is currently in the process of 
being converted to an electronic database and GIS format. Site and inventory data for Elko 
County was previously entered into the statewide GIS, and into a Microsoft Access® database. 
The database contains fields and codes identical to the IMACS site record. Spatial and 
database information for sites and inventories lying within the White Pine county and 
Lincoln county portions of the study area were compiled as part of this project. 


Several steps were involved in data compilation for the study area within Lincoln and White 
Pine counties. First, archival USGS maps (7.5 and 15 minute quadrangles) containing site 
and inventory locations were scanned at the UNLV Harry Reid Center archive. Those quads 
were then geo-referenced to UTM Zone 11; NAD 27 coordinates. Each site and inventory 
marked on the maps was digitized. Any sites smaller than 2.5 acres in extent were digitized 
as point features using GIS software; linear sites were digitized as lines; all other sites were 
represented as polygons. Similar digitizing rules were applied to inventoried areas. Site and 
inventory metadata consisting of map source, entry dates and accuracy or error flags were 
appended to attribute tables for each shape. 


Site data from records predating IMACS (1982) proved to be somewhat inconsistent. 
Likewise, early investigations are generally less complete than more recent ones and the 
survey methods used at the time varied considerably. To control for variability in survey 
method and site reporting, assemblage and administrative site data were entered only for 
those sites occurring within inventories with a cumulative extent greater than 640 acres. Size 
criteria assured relatively uniform reconnaissance and reporting technique and constrained 
site vs. non-site analysis of the landscape within consistent parameters. 


Sites were selected by intersecting inventory area with site location. Site records were 
assembled from archives at the BLM Ely Field Office, UNLV Harry Reid Center and the 
Nevada State Museum. Administrative and assemblage data were compiled in an Microsoft 
Access® database using the IMACS encoding format, then linked to the spatial data in the 
GIS attribute tables. Like the Utah data, shapefiles were transformed into a single polygon 
layer by buffering points and lines into a synthetic polygon shape, then merging those with 
the existing polygon shapefile for analytical purposes. Site centerpoints were also calculated 
for each feature for use if point analysis was required. 
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' Study Area Data Files 


After site data from all three states were assembled, GIS shapefiles were merged into a single 
analytical theme and joined to respective site assemblage data. Since assemblage data was 
reported in slightly different format for each state, attribute fields were reformatted to 
indicate presence or absence of specific artifact types or general classes, feature types, and 
temporal affiliation. The resulting table produced comparable data attributes for all site 
records. It was used to identify historic and prehistoric site affinity and created a baseline for 
archaeological and anthropological site analysis. Figure 4.1 depicts the distribution of site 
center points across the study area. Inventories greater than 640 acres are shown within the 
Nevada and Utah data set. 


Evidential Themes 


Landscape level analysis required the compilation of a number of environmental data sets or 
evidential themes that could be used with the site data to construct a probability model. Data 
sets compiled for the project area included slope, vegetation, landform, and hydrology. A 
roads layer was compiled for historic resource analysis. GIS layers pertaining to potential 
marsh habitat were also derived as a means to address research questions relating to 
prehistoric land use. 


Slope 


Slope was derived from the USGS National Elevational Data set (NED). The 30 meter NED 
was clipped to each analytical unit within the project area and slope was calculated for each 
cell, and then converted to a slope grid. For analytical purposes, slope was divided into five 
classes: 0-5 degrees, 5-15 degrees, 15-30 degrees, 30-45 degrees and greater than 45 degrees. 
The NED was also used to create shaded relief maps for use as background graphic in each 
of the analytic units. 


Vegetation 


Vegetation layers were derived from Fire Sciences Laboratory, Rocky Mountain Research 
Station, Potential Natural Vegetation Groups (Schmidt et al. 2002). This is coarse-scale data 
that were developed as part of a national level, fire-planning model. Vegetation data was 
refined to match terrain using a 500 meter Digital Elevation Model, 4" Code Hydrological 
Units and Ecological Sub-regions (Bailey’s Sections). Classifications follow Kiichler (1975) 
descriptions for ECO Region 4 (Table 4.2). 


Landform 


In order to derive a general characterization of landform within each analytic unit, the NED 
data set was reclassified into three ranges of slope that roughly approximate flats, piedmont 
and mountainous areas. Flats comprise all slopes between 0 and 3%; piedmont lies between 4 
and 10%; and mountains are all slopes above 10%. The resulting classes approximate 
elevational rings of valley bottom, alluvial fan and upland slopes for each analytic unit. 
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Table 4.2 
Potential Natural Vegetation Groups ECO Region 4 (after Kichler; 1975) 
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Springs and Streams 


A hydrologic layer consisting of springs and streams was compiled for each of the analytic 
units. Source data was derived from USGS 1:100,000 Digital Line Graphs (DLG) clipped to 
the project area then buffered at intervals of 200, 400, 1000 and 2000 meters. Buffered 
shapes were then converted into grids for each analytic unit. Both intermittent and perennial 
stream classes are included in the data set, since present intermittent water courses may have 
been more productive prehistorically. 


Potential Wetlands 


The extent of potential marsh habitat was derived from the U.S. Soil Conservation Service 
STATSGO State Soil Geographic database. The STATSGO database was designed for use as 
a regional, multi-state resource planning, management and monitoring tool. Soil data is 
derived from generalized information provided in the county-wide soils database and 
extrapolated to 1:250,000 scale USGS quadrangles. STATSGO data sets include fields 
relating to soil class, structure, texture, engineering capabilities, suitability for agriculture, 
and potential for various rangeland habitat types. STATSGO databases were queried for soils 
with the potential to sustain wetland plants and the potential to sustain wetland wildlife. The 
results were used as a proxy for potential wetlands. Those shapes were then buffered at 1000, 
3000 and 5000 meter intervals for analytical purposes and then converted to grids. 


Roads 


The roads layer was extracted from U.S. Census Bureau 2000 Tiger/line files. Data was 
derived from a generalized 1:100,000 base layer. Line data was then buffered to 200, 400, 
and 1000 meter widths for analytic purposes. 


Analytic Methods 


Cultural resource inventory data sets allowed for multiple approaches be used to construct a 
management model. Inventoried areas provided a controlled setting where both site and non- 
site data can be assessed. Within the site/non-site parameters chi-square analysis could also 
be conducted to validate predictive patterns observed in the calculated weights tables. 


Weights tables were compiled in Spatial Data Modeler using sites within inventoried areas 
as a training point theme and inventory extent as a mask over all evidential themes, Unit area 
settings suggested by Spatial Data Modeler vary according to analytic unit size. The 
suggested unit area compensates for variation between study area cell size and output cell 
size of the evidential themes (Suggested Value= (total Study Area / total Training Points) / 
40). Default settings for most of the analytic units ranged from 0.20 to 0.30 square kilometers 
(447.2 or 547.7 meter grid). To maintain consistency within each analytic unit, the unit area 
was arbitrarily set to 500 meter cells (0.25 square kilometers). Multiple training points within 
a cell greatly inflate prior probabilities since probability is evaluated as a deviation from the 
normal distribution of one training point per unit area. SDM will automatically weed or 
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remove any duplicate training points within a cell so that there are no more than one training 
point (site) per unit area. (< 


Once unit area and training point parameters are set, Spatial Data Modeler calculates a 
weight table for each evidential theme. The resulting contrasts (weight * - weight’) indicate 
the relative strength of each predictive class. 


To test the efficacy of the weights calculations and to aid in the selection of predictive classes 
for creation of a final response theme, a chi-square test was run with the inventoried site data 
set against the evidential themes. To create a site/non-site matrix, the project area was 
arbitrarily gridded into 250 meter square cells and a centerpoint was calculated for each cell. 
Centerpoints were clipped to the analytic unit, then 1 apain clipped so that only grid points 
within inventoried areas remained. Using ArcView” Spatial Analyst, any grid point within 
100 meters of a site polygon was selected and saved as a site training point. The selected 
subset was switched, and all remaining grid points were saved as a non-site theme (Figure 


4.2). 


With the Spatial Data Modeler area unit set to 250 meter cell size, weights were calculated 

using both site and non-site training point themes. Resulting contrasts were compared with 

the previous run of weeded, inventoried sites. Classes with the highest contrasts in both the 

250 meter grid site and weeded site weights tables were chosen for validation using the chi- 

square test. Evidential class with the highest contrast was tabulated against site and non-site 
occurrences (Table 4.3). A chi-square above 3.84 was considered significant at 1df. Lf chi- 

square testing confirmed the contrast as predictive, that class was chosen as “inside” the “ 
pattern. 


Response themes, using all sites (weeded) within each analytic unit, and the predictive 
classes were then run for each of the hydrographic units that contained inventory themes. The 
normalized posterior probability was then reclassified to reflect high, moderate and low 
probability of site occurrence. Summary tables of sites within each probability zone were 
compared with the results from areas of previous inventory. The probability model was 
considered accurate if highest site frequencies were associated with areas of high and 
moderate probability. Since the Idaho data lacks spatial data for inventories, the comparisons 
allowed us to assess the feasibility of using site center points regardless of inventory status as 
valid training points for pattern prediction. 
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Table 4.3 
Chi-Square Analytic Format 


Inside Class 
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Vv. ANALYSIS 


Of the 83,000 square kilometers within the study area, 78,747 square kilometers were 
evaluated as part of the probability model; 4350 square kilometers of land in the Upper 
Snake analytic unit not under Bureau of Land Management control were excluded. (Table 
5.1) Systematic inventories have been conducted over approximately 4% of the model, and 
inventories greater than 640 acres in extent comprise 80% of that area. A total of 5284 sites 
are reported within the model area, 1819 of them fall within the larger inventory blocks. The 
following chapter describes respective hydrologic units and presents results of the probability 
model for each analytic unit within the model area. 


PILOT/THOUSAND SPRINGS VALLEY ANALYTIC UNIT 
Analytic Unit Description 


The Pilot Springs/Thousand Springs analytic unit covers approximately 1.1 million acres 
(1785 mi’)/.4 million hectares (4623 km/), It lies in within the northeastern comer of Nevada 
with a small portion falling within western Utah. The analytic unit lies within the Great Basin 
region, but with drainage eastward into the Great Salt Lake Desert and Bonneville Basin is 
considered a sub-unit of the Great Salt Lake hydrographic unit. The upland characterization 
of this analytic unit drove the decision to analyze it separately from the larger Great Salt 
Lake analytical unit. (Figure 5.1) 


Several small valleys and basins comprise the Pilot/Thousand Springs analytic unit. 
Thousand Springs Valley and Pilot Springs Valley are the most predominate, covering a 
major portion of the analytic unit. Toano Draw slopes northward into Thousand Springs 
Valley and Tecoma Valley extends north from Pilot Springs Valley. A number of relatively 
low ranges and mountains define the limits and interior of the analytic unit. The Toano 
Range, Pequop Mountains, Windemere Hills and the Snake Mountains define the 
southwestern extent of the Pilot/Thousand Springs area. Knoll Mountam and Cedar Mountain 
mark the hydrographic units northern extent, while the Delano Mountains, Pilot Range, and 
Leppy Hills form an eastern boundary. Ninemile Mountain and Murdock Mountain separate 
Toano Draw and Thousand Springs Valley from the eastern valleys. Elevations of the 
surrounding mountains are relatively low, extending between 2200 and 2700 meters ams. 


All valleys within the Pilot/Thousand Springs analytic unit are externally drained. Thousand 
Springs Creek Flows north and eastward from Toano Draw and the Snake Range around 
Ninemile Mountain, then southeasterly through Tecoma Valley into the northwestern uplands 
of the Great Salt Lake Desert. Pilot Springs Creek drains southward through Pilot Springs 
Valley then terminates in an extensive sand sheet and dry flat between the southern extent of 
the Pilot Range and the Leppy Hills. The southern extent of Pilot Springs Valley lies at 1340 
meters amsl, just above the Gilbert Shoreline of Lake Bonneville. Toano Draw and Thousand 
Springs Valley lie at elevations between 1800 and 1600 meters. As it drains through Tecoma 
Valley, Thousand Springs Creek attains an elevation of 1400 meters amsl, 





Table 5.1 
Summary Sites and Inventories Within Each Analytic Unit 


Model Area [Total# inventory |__| # 640+ Acre 
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Figure 5.1 Pilot/Thousand Springs Valley Analytic Unit 
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Vegetation within Pilot/Thousand Springs analytic unit is primarily sagebrush with juniper 
and juniper/pinyon forest on mountain slopes. Barren areas occur in the southern dunes and 
flats while desert shrub communities in the lower portions of Pilot Springs Valley and 
Tecoma Valley. 


Analytic Results 
Prehistoric Evidential Themes 


Of the 4622 square kilometers within the Pilot/Thousand Springs Valley analytic unit 
approximately 3.5% (164 km”) have been assessed by inventories larger than 640 acres. Four 
hundred sixty prehistoric sites are reported within those inventories, while 697 sites are 
reported within the analytic unit as a whole. (Table 5.2) (Figure 5.2) 


Sampling within each of the evidential classes is relatively consistent. The juniper steppe 
vegetation zone is less than 1 square kilometer in extent and has not been sampled. Less than 
2% of the juniper/pinyon zone has been inventoried. Areas lying more than 1000 meters from 
streams and between 3000 and 5000 meters from potential wetlands have also been poorly 
sampled, 


Calculated weights for each evidential theme suggests that a predictive pattern for sites 
occurs within the desert shrub vegetative community, within 1000 meters of potential 
wetlands and within piedmont slopes.(Table 5.3) Positive contrasts for slope and distance to 
springs or streams are inconsistent across analytic runs, and calculated chi-squares suggest a 
normal distribution of sites within high contrast classes. 


Within vegetation evidential themes, desert shrub is the only class with a high contrast. 
While 10% of the sites lie within the juniper/pinyon zone, distribution of sites is less than 
anticipated for a positive pattern association. (Figure 5.3) 


Potential wetland areas within the Pilot/Thousand Springs analytic unit are relatively few, as 
reflected by the cumulative extent of those areas lying outside of the 5000 meters buffered 
area. Areal extents of the three buffered zones are approximately equal, and contrast is 
uniformly high and strongly predictive for the 0-1000 meter buffer. (Figure 5.4) 


When only inventoried areas are considered, the piedmont is the most predictive class for 
sites. An analytic run using all sites identifies flats as the most predictive class, but by 
controlling for inventoried space, the number of sites within that area is reduced by almost 
42%. By contrast, 78% of all sites within the piedmont landform are accounted for by 
inventories greater than 640 acres in extent. (Figure 5.5) 


- Prehistoric Predictive Response 
Posterior probabilities generated within the response theme for the Pilot/Thousand Springs 


Valley analytic unit cluster within three groups with breaks at 0.072 and at 0.044. The prior 
probability for a normal distributional pattern is 0.029, well below the (Table 5.4) (Figure 
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Table 5.2 
Pilot/Thousand Springs Valley Analytic Unit Inventory Summary 
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Figure 5.3 Pilot/Thousand Springs Valley Analytic Unit Predictive Pattern - Vegetation 
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Figure 5.4 Pilot/Thousand Springs Valley Analytic Unit Predictive Pattern - Potential Wetland 
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Figure 5.5 Pilot/Thousand Springs Valley Analytic Unit Predictive Pattern - Landform 
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Pilot/Thousand Springs Valley Analytic Unit Initial Response 
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5.6) discernable break for low probability. With breaks at those points, well over one-half of 
the sites within the analytic unit falls within the area of lowest probability. In an attempt to 
capture additional sites and balance the distribution of sites within each probability zone, 
breaks were re-drawn at 0.044 and 0.029, just above the intersection with prior probability. 
(Table 5.5) (Figure 5.7) Results of the model derived from those parameters were tallied and 
show an almost even distribution of sites within high to moderate and low probability areas. 
(Table 5.6) The results appear to be biased by a relatively high frequency of sites (training 
points) lying within the large area of flats and sagebrush, with a corresponding low weight 
and contrast relating to a normalized distribution. (Figure 5.8) 


The response table presents a tally of presence or absence of predictive evidential classes, 
then recalculates probabilities based upon the area and number of training points within each 
row of tabulated intersections. If large areas contain a proportional number of sites, 
probabilities will by definition remain near the prior probability. Likewise, negative weights 
and negative contrasts will still retain their lower probabilities in the response theme. 


The logic behind the response theme is that as intersecting predictive themes overlap, 
corresponding probabilities validate the predictive relationship within each defined class. 
Probability based correlations fail when a significant number of the evidential classes exhibit 
negative contrasts as a result of lower than expected frequencies within disproportionately 
large areas. 


In order to derive a version of the response theme based upon the overlap of predictive 
evidential classes, evidential themes were reclassified using the binary values assigned to 
inside or outside pattern within each theme. Using Spatial Analyst®, a new class was 
calculated by combining each of the predictive layers into a single class. Rows containing 1, 
for presence within the predictive pattern, were totaled, with results ranging from 0, no 
overlap present to 3, all three themes intersect. Those results were then re-classified into low 
(0 overlap), medium (1 class present) and high (2 or 3 classes present) sensitivity zones. 


The resulting response presents a better fit of probability layers to the actual site area. Total 
area varies between the two different response runs due to grid variation within the 
vegetation evidential theme. (Table 5.7) (Figure 5.9) The distribution of sites within high 
and medium probability zones comprises more than 70% of the total site area within 55% of 
the total model area. Ratios of site area to model and inventory area exhibit the same trend, 
with highest ratios descending significantly from high to low sensitivity zones. 


Historic Evidential Themes 


Sixty-nine historic sites are recorded within the Pilot/Thousand Springs Valley analytic unit. 
Of those, only 19 (28%) fall within inventories greater than 640 acres in extent. (Table 5.2) 
(Figure 5.10). Weights tables for the historic evidential themes, indicate varying positive 
contrasts within buffered distances to roads and water. Chi-square for roads is significant at 
the 400 meter buffer of inventoried sites, but is not significant for distance to water. When 
buffer areas for roads between 0 and 400 meters are combined, chi-square remains 
significant. 





Figure §.6 PilotThousand Springs Valley Response Breaks 
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Table 5.5 
Pilot/Thousand Springs Valley Analytic Unit Observed Response 
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Figure 5.7 Pilot/Thousand Springs Valley Observed Response Breaks 
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Table 5.6 
Pilot/Thousand Springs Valley Analytic Unit Model Summary Prehistoric Observed Response 
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Figure 5.8 Pilot/Thousand Springs Valley Analytic Unit Observed Probability - Prehistoric 
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Table 5.7 
Pilot/Thousand Springs Valley Analytic Unit Model Summary Prehistoric Composite 
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Note: Total area may vary between response and composite analysis due to grid variation within the vegetation 
evidential theme. 





Figure 5.9 Pilot/Thousand Springs Valley Analytic Unit Composite Probability - Prehistoric 


[_] Pilot Analytic Unit 
+ Prehistoric Sites (inventoried) 
o Prehistoric Sites 


Probability 10 0 10 20 Miles N 
[_] Low | A 
Wi Medium 

High 20 0 20 40 Kilometers 


——e 


| | No Data 





> 





Fy a € 
‘ Bey % i 
y , wed 


Figure 5.10 Pilot/Thousand Springs Analytic Unit - Inventories and Historic Sites 
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Historic Predictive Response 


Since a response theme cannot be created with only one class inside the pattern, the 0-400 
meter buffer from water was used as a predictive theme. The resulting grid classifies 
probability as high or low; area within 400 meters is high, greater that 400 meters is low. 
(Figure 5,11) Summary tables reflect the expected site distribution with 64% of the analytic 
unit comprising the low probability zone, while slightly more than 8% of all sites fall within 
that area. More than 90% of all sites and 85% of inventoried sites fall within the high 
probability zone. (Table 5.8) (Figure 5,12) 


RUBY/LONG VALLEY ANALYTIC UNIT 
Analytic Unit Description 


The Ruby/Long Valley analytic unit is the hydrographic unit within the GBRI study area. It 
shares its eastern boundary with the Spring/Steptoe Valley Analytic unit, and its northem 
extent with the Pilot/Thousand Springs Valley analytic unit (Figure 5.13). In addition to 
Ruby and Long valleys, the analytic unit includes Clover Valley and Independence Valley in 
the north along with Butte Valley and Jakes Valley in southeast. The analytic unit covers 
approximately 2.6 million acres (4095 mi’)/1.0 million hectares (1060 km”). Bounding 
ranges of the hydrographic unit include the White Pine Range, Ruby Mountains, and 
Humboldt Range to the west, Wood Hills and Windemere Hills to the north and the Pequop 
Mountains, Cherry Creek Range and Egan Mountains to the east. The Maverick Springs 
Range, Butte Mountains and Medicine Range provide a barrier between Ruby/Long Valley 
in the western portion of the analytic unit and Butte/Jakes Valley in the east. 


Elevations of the Ruby Mountains and Humboldt Range exceed 3000 meters amsl. Northern 
ranges are lower, averaging 2700 meters ams] while southern bounding ranges and interior 
ranges extend to 2800 meters amsl. Likewise, valley floor are relatively high averaging 2000 
meters in elevation with valley floors between 1850 and 1800 meters. 


Hydrologically, each of the valleys within the analytic unit is internally drained. The Franklin 
River and the Ruby Marshes, consisting of Ruby Lake and Franklin Lake are the major 
hydrographic features within Ruby Valley. Snow Water Lake serves as a major hydrologic 
collection point for Clover Valley. Bounding mountains of the remaining valleys provide 
ample perennial flow, but all terminate in dry flats at the valley bottom. Faulting has 
produced numerous springs along the steeper eastern escarpment of the bounding and interior 
mountain ranges. 


Vegetation is similar to that in Spring/Steptoe Valley. Limber pine and alpine vegetation 
occurs on the highest slopes, with juniper/pinyon woodlands on lower more protected slopes. 
Riparian meadows and wetland habitat dominates the area of perennial lakes and marshes, 
while sagebrush is the dominant vegetation on the piedmont and upper valley slopes. Lowest 
portions of the valley floor consist of desert shrub communities while dry flats and valley 
bottomland is sparsely vegetated. 
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Figure 5.11 Pilot/Thousand Springs Valley Analytic Unit Predictive Pattern - Roads 
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Table 5.8 
Pilot/Thousand Springs Valley Analytic Unit Model Summary Historic Composite 


























TC High (1) (0-400) | Medium | Low (0) (>400 
Model area (m’) 1652891008.00] __—=—s'|_—_—-2969638912.00] 4622529920.00 
Modelarea i?) | ___tesz.eal 0.00] eeneal_ 422.59 
% Model area 35.76% 64.24% 100.00% 
All sites area (m*) as 2138905,0000]  ==s| Ss: 192888.7188] ——«2331793.72 
All sites area (km?) pe: 51 HT aS 
91.73%] 0.00% 8.27% 400.00% 
All site area/ model area__ 0,0013 
inventory area (m’) 6485455200] Ss | =~: 99525368.00] _ 164379920.00 
Inventory area (km’) Bae ST" ee eae 164.38 
39.45%] __0.00%| 60.55% 100.00% 
% Model area inventoried 392%| 0.00%| 8.35%] (3.56% 
inventory sites area (m”) 274065.84, 49153.11 323218.95 
inventory sites area (em?) | ——oz7| ooo] _—o.ast—— 0.22 


0.27 
% Inventory site area_ 0.009 15.21%| __-100.00% 


Inv site area / inv area 0.0042] 0.0000 0.0005 





Note: Total area may vary between response and composite analysis due to grid variation within the vegetation evidential theme. 





Figure 5.12 Pilot/Thousand Springs Valley Analytic Unit Composite Probability - Historic 
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Analytic Results 
Prehistoric Evidential Themes 


Of 10,606 square kilometers in the Ruby/Long Valley analytic unit, approximately 927 
square kilometers, 8.5% of the total area, have been inventoried. (Figure 5.14) Six hundred 
thirty-eight sites are reported as part of inventories greater than 640 acres, 973 sites are 
identified within the entire analytic unit. (Table 5.9) 


All analytic classes, were sampled during previous inventories. The Ruby Marshes and 
surrounding marsh and wetland habitat have been extensively investigated, and unlike most 
of the analytic units, steeper slopes have been more intensively examined. 


Weights of evidence tables identify classes within each evidential theme that lie “inside” the 
predictive pattern. (Table 5.10) Normalized contrast for meadows are highest in all runs of 
prehistoric sites. Other vegetation classes display negative or very low positive contrasts. 
Sagebrush and water have relatively high contrasts when all sites are considered. Meadows 
within the analytic unit cover less than 2% of the entire area and are considered part of the 
marsh environment. (Figure 5.15) 


Contrasts for distance to springs and streams are consistently high for inventoried areas 
between 1000 and 2000 meters from that class of water. Buffered areas from 200 to 1000 
meters from a water course are consistently sampled, but reveal lower than expected or 
marginal site frequencies. The lowest contrasts are evident at distances more than 2000 
meters from springs and streams. (Figure 5.16) 


Proximity to wetlands is highly predictive within the Ruby/Long Valley analytic unit. 
Highest contrasts are evident within 1000 meters of potential wetland habitat, while areas 
lying more than 3000 meters from that zone show a negative correlation with normal site 
distribution. Proximity to wetlands correlates well with vegetation contrasts. (F igure 5.17) 


Slopes between 0 and 5 degrees are highly predictive for sites within this analytic unit. 
Nearly two-thirds of the inventoried area occurs on flat slopes, and positive contrasts are 
evident on slopes up to 15 degrees. Slopes above 15 degrees uniformly exhibit a negative 
contrast. (Figure 5.18) 


Landform strengthens the relationship of slope as a predictive theme in the Ruby/Long 
Valley analytic unit. When all sites are considered, both flats and piedmont have a high 
predictive contrast, while inventoried areas show highest contrasts within the piedmont. Chi- 
square statistics confirm a non-random distribution of sites on the piedmont. (Figure 5.19) 


Prehistoric Predictive Response 
Normalized posterior probabilities were used as a means to evaluate tabular results from the 


response theme generated for the Ruby/Long Valley analytic unit. (Figure 5.20) Prior 
probability for the response theme was set at 0.0181 and observed breaks within normalized 
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Figure 5.15 Ruby/Long Valley Analytic Unit Predictive Pattern - Vegetation 
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Figure 5.16 Ruby/Long Valley Analytic Unit Predictive Pattern - Streams and Springs 
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Figure 5.17 Ruby/Long Valley Analytic Unit Predictive Patter - Potential Wetland 
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Figure 5.18 Ruby/Long Valley Analytic Unit Predictive Patter - Slope 
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Figure 5.19 Ruby/Long Valley Analytic Unit Predictive Pattern - Landform 
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Figure 5.20 Ruby/Long Valley Analytic Unit Observed Probability - Prehistoric 
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posterior probability were set at 0.182 and 0.0033. (Table 5.11) (Figure 5,21) Highest 
probabilities for encountering sites occur when evidential classes identified as inside the 
predictive pattern intersect. Combinations of proximity to wetlands, springs and streams, 
vegetation, and slope have the highest posterior probabilities, but combinations of three or 
more evidential themes are also common within the range of moderate probabilities and two 
or more combinations of evidential themes within the low probability area. 


Summary tables for the Ruby/Long Valley analytic unit response theme shows that while 
33% of the area is classified as low probability, over 57% of the sites occur within that zone. 
A similar pattern occurs for inventoried sites. Table 5.12 The higher frequency of training 
points within the low probability area creates a normal distribution and with it, a posterior 
probability lower than the prior distribution. 


To realign response themes so that they better correspond with predictive patterns within 
evidential themes, new probability areas were calculated by totaling binary theme values of 
each theme. Since all five themes contained predictive classes, additive scores ranged from 0 
to 5. Probability classes were grouped into three classes, 0-1 low, 2 medium, and 3-5 high. 
(Table 5.13) Site area distributions within newly defined probability zones provide a better 
fit of the data. Site densities are highest in the high probability zones and lowest in low 
probability areas. Slightly more than 11% of the site area occurs within the low probability 
zone, which accounts for slightly less than 20% of the analytic unit. (Figure 5.22) 


Historic Evidential Themes 


One hundred fifty-seven historic sites are reported within the Ruby/Long Valley analytic unit 
and 81 of these are located within 640 acre or larger inventory units, (Table 5.9) (Figure 
5.23) The area within 200 meters of roads, and within 200 meters of streams or springs, 
revealed the highest contrast within historic evidential themes. {Figure 5.24) (Figure 5.25) 
Nearly 75% of the inventoried sites (61) lie within 200 meters of roads, while 32% of 
inventoried sites (26) lie within 200 meters of potential water sources. Distances greater than 
200 meters are uniformly less predictive. (Table 5.14) 


Historic Predictive Response 


Historic response themes generated for the Ruby/Long Valley analytic unit show three 
possible breaks in the posterior probabilities; 0.014 to 0.008, 0.008 to 0.003 and 0.003 to 
0.0009, with a prior probability set at 0.003. (Table 5.15) (Figure 5.26) Lower relative 
contrasts for proximity of sites to potential water sources create a cluster of training points 
with posterior probabilities below the prior expected value. Summary tables show that the 
resulting probability map meets expectations for site density in the low probability area. 
Seventy-five percent of the analytic unit comprises the low sensitivity zone, and 25% of the 
sites fall within this area. (Table 5.16), Fourteen training points associated with proximity to 
water are associated with the low probability zone. The medium probability zone is relatively 
small and contains a single set of training points associated only with proximity to roads. 
(Figure 5.27) 
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| /Long Valley Analytic Unit Prehistoric Response 
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Figure 5.21 Ruby!Long Valley Analytic Unit Response 
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Table 5.12 
Ruby/Long Valley Analytic Unit Model Summary Prehistoric Response 
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Table 5.13 
Ruby/Long Valley Analytic Unit Model Summary Prehistoric Composite © 
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Figure 5.22 Ruby/Long Valley Analytic Unit Composite Probability - Prehistoric 
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Figure 5.23 Ruby/Long Valley Analytic Unit - Inventories and Historic Sites 
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Figure 5.24 Ruby/Long Valley Analytic Unit Predictive Pattern - Roads 
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Figure 5.25 Ruby/Long Valley Analytic Unit Predictive Pattern - Water 
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Table 5.14 
Ruby/Long Valley Analytic Unit Historic Evidential Theme Weights/Chi-Square 
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Table 5,15 
Ruby\Long Valley Analytic Unit Historic Response 
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Figure 5.26 Ruby/Long Valley Analytic Unit Historic Response 
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Table 5,16 
Ruby/Long Valley Analytic Unit Model Summary Historic Response 
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Ruby/Long Valley Analytic Unit Model Summary Historic Composite 
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Figure 5.27 Ruby/Long Valley Analytic Unit Observed Probability - Historic 
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Using similar methods employed to clarify the prehistoric predictive response, the historic 
response theme was re-classified using only the intersection of predictive evidential classes. 
A more balanced distribution is evident within the reclassification. (Table 5.1 6) Forty three 
percent of the analytic unit is classified as high or medium probability, with 84% of the 
inventoried site area falling into that classification. Sites falling within 200 meters of water 
previously associated with the low probability zone are now included within the medium 
classification. (Figure 5.28) 


SPRING/STEPTOE VALLEY ANALYTIC UNIT 
Analytic Unit Description 


The Spring/Steptoe Valley analytic unit lies to the east of the Ruby/Long Valley analytic unit 
and includes Spring Valley and Steptoe Valley to the south, and Goshute Valley and 
Antelope Valley to the north. The analytic unit covers approximately 3.4 million acres (5323 
mi’)/1,3 million hectares (13787 km’). (Figure 5.29) Topography is typical of north-trending 
grabens within the Great Basin. High bounding ranges create an orthographic effect on 
precipitation patterns depositing more moisture along west-facing slopes. Steeply faulted 


bounding ranges produce numerous springs along eastern pediment slopes. 


Steptoe Creek and Duck Creek are the major hydrologic features along the western side of 
the analytic unit. Both drain northward into Goshute Lake at the southern end of Goshute 
Valley. Spring Valley Creek is the major drainage in the eastern portion of the analytic unit. 
It flows north through Spring Valley, then terminates in a large depression and dune field 
south of Spring Creek Flat. Antelope Valley is relatively dry. Drainages flow from the 
surrounding mountains to the valley floor. Numerous spring complexes occur within Spring 
Valley, especially along the toe of western piedmont slopes. Marshes and ponds are present 
in Steptoe Valley along Steptoe Creek southeast of Ely and west of McGill. 


Elevations of the valley floors within the Spring/Steptoe analytic unit are relatively high, 
ranging from 1900 meters in the south to 1750 meters in the north. The Pequop Mountains 
and Toano Range bound the hydrographic unit in the north, Cherry Creek and Egan Range on 
the west and the Snake Range and Ferber Hills to the east. The Schell Creek Range separates 
Steptoe and Spring Valleys. Wheeler Peak (3952 meters) in the Snake Range is the highest 
peak within the analytic unit. Mountain elevations are highest in the southern portion of the 
analytic unit, averaging 3500 meters. Norther ranges average approximately 2500 meters in 
elevation. 


Vegetation is typical of the Great Basin. Highest elevations are dominated by alpine 
vegetation including limber and bristlecone pine; juniper/pinyon forest covers more 
temperate lower slopes. The sagebrush zone dominates open pediment slopes and is replaced 
by desert shrub communities on the lower flats. Depressions and valley bottoms are sparsely 
vegetated, 
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Figure 5.28 Ruby/Long Valley Analytic Unit Composite Probability - Historic 
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Figure 5.29 Spring/Steptoe Valley Analytic Unit 
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Analytic Results 
Prehistoric Evidential Themes 


Approximately 387 square kilometers (2.8%) of 13,787 square kilometers within the entire 
Spring/Steptoe Valley analytic unit were inventoried. (Figure 5.30) Eight hundred twenty 
three sites were identified within the analytic unit; 410 (50%) of those were recorded in 
inventories larger than 640 acres in extent. Table 5.17 


Larger surveys have occurred within a sample of all analytic classes within the analytic unit. 
However, vegetation zones within Great Basin pine, and barren areas along with mountain 
landforms and slopes greater than 15 degrees are under sampled. Great Basin pine and barren 
vegetation zones respectively comprises 0.9% and 1.2% of the of the entire model area. 
Approximately 28% of the analytic unit comprises the mountain landform. Thirty-nine 
percent of the analytic unit consists of slope greater than 15 degrees. 


Evidential theme classes with highest contrasts and correspondingly significant chi-square 
results were identified as lying “inside” the predictive pattern. Table 5.18 Highest contrast 
for vegetation within inventoried sites was associated with sagebrush, and the associated chi- 
square was also considered significant. Sagebrush retains its high contrast when all sites are 
considered. A high negative contrast within the juniper/pinyon class along with a high 
positive contrast for non-sites within that class strongly suggests a lower than expected 
relationship between sites and the juniper/pinyon zone. (Figure 5.31) 


Highest contrast for distance to water is within the 1000 to 2000 meter buffer band. While 
greater numbers of sites are located between 0 and 1000 meters of streams and springs, their 
numbers reflect a normal distribution in relation to area. The band between 200 and 400 
meters shows the strongest negative contrast with fewer sites than would normally be 
expected within that region. (Figure 5.32) 


Very few potential wetland environments are present within the Spring/Steptoe Valley 
analytic unit and less than 20% of the analytic unit lies within 5000 meters of this evidential 
theme. Not surprisingly, highest site contrasts within inventoried areas are within the 
buffered class that lies more than 5000 meters from potential wetlands. When all sites are 
considered, areas more than 5000 meters from wetlands have a negative contrast, indicating a 
weak correlation with sites. Since the results of the weights calculations are inconclusive, no 
classes within potential wetlands were selected as most predictive for analysis. 


Weight calculations for slope in all runs identifies gradient between 0 and 5 degrees as 
having the highest contrast. Likewise, chi-square for slopes between 0 and 5 degrees meets 
the critical value for non-random distribution. While sampling discrepancy may marginally 
effect contrast values, the high frequency of sites both within inventoried samples and when 
all sites are considered suggests that slopes of 0 to 5 degrees are “inside” the predictive 
pattern. (Figure 5.33) 
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Figure 5.30 Spring/Steptoe Valley Analytic Unit - Inventories and Preh 
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Table 5.17 
ns Valley Analytic Unit Inventory Summary 
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Figure 5.31 Spring/Steptoe Valley Analytic Unit Predictive Pattem - Vegetation 
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Spring/Steptoe Valley Analytic Unit Prehistoric Evidential Theme Weights/Chi-Square 


Table 5.18 
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Figure 5.32 Spring/Steptoe Valley Analytic Unit Predictive Pattern - Streams and Springs 
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Figure 5.33 Spring/Steptoe Valley Analytic Unit Predictive Pattern - Slope 
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Flats exhibit the highest contrast for landforms in all analytical runs. Frequency of sites is 
highest in this class and suggests a greater than normal distribution of sites. In contrast, the 
negative contrast for piedmont with adequate sampling suggests that fewer sites than 
expected occur within that class. (Figure 5.34) 


Prehistoric Predictive Response 


Normalized posterior probabilities were used as a means to evaluate tabular results from the 
response theme generated for the Spring/Steptoe Valley analytic Unit. (Figure 5.35) Prior 
probability for the response theme was set at 0.0122 and observed breaks within normalized 
posterior probability were set at 0.016 and 0.0095. (Table 5.19) (Figure 5.36) Highest 
probabilities for encountering sites occur as evidential classes identified as inside the 
predictive pattern intersect. Three or more combinations of flats and 0 to 5 degree slope, with 
sagebrush or 1000 to 2000 meters from water have the highest probability scores and the 
highest frequency of training points. As intersecting conditions decrease, probability 
becomes moderate and where only a single evidential condition is met, probabilities fall 
below the prior value of 0.012 and are low. Analysis of the response themes show that by 
area, the highest proportion of sites fall within high and medium probability zones. (Table 
5.20) High probability zones have been more intensively surveyed than zones of medium or 
low probability and the lowest proportion of inventoried site by area occurs within the 
medium probability zone, When high and medium probability zones are combined, however, 
over 80% of the inventoried site area falls within that area, About 66% of all site areas fall 
within the high and moderate probability zones. 


When probabilities are recalculated by composite predictive class, the extent of high and 
medium probability areas are increased and of low probability area is decreased. (Table 5.20) 
Twenty-nine percent of the analytic unit lies within the high probability zone and 49% of all 
sites fall within that area. Forty-four percent of the analytic unit is classified as medium 
probability and 44% of the sites occur there, while 19% of all site areas fall within the 
remaining 27% of the analytic unit. Within inventoried site areas, almost 73% of the sites 
fall within high probability areas, 19% in medium probability zones and 8% in the low 
probability zone. (figure 5.37) Over 56% of the inventories have been conducted within 
areas identified as high probability and most of the site area lies within that zone. 
Correlations between site density and high probability areas may be biased by sampling 
within the Spring/Steptoe Valley analytic unit 


Historic Evidential Themes 


One hundred forty-four historic sites are included in the analysis of the Spring/Steptoe Valley 
analytic unit. Seventy-nine of these have been identified within inventories greater than 640 
acres in extent. (Table 5.17) (Figure 5.38) 


Contrasts within buffered classes of roads and water sources are easily discernable. Highest 
contrasts are evident within 200 meters of existing roads, and weights for non-sites exhibit a 
negative contrast within the same buffer. (Table 5.21) ( Figure 5.39) Other contrasts for 
buffered distances to roads are either negative or lightly positive. Highest contrast and 








Figure 5.34 Spring/Steptoe Valley Analytic Unit Predictive Pattern - Landform 
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Figure 5.35 Spring/Steptoe Valley Analytic Unit Observed Probability - Prehistoric 
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Table 5.19 | 
& Spring/Steptoe Valley Analytic Unit Prehistoric Response 
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Table 5.20 
Spring/Steptoe Analytic Unit Model Summary Prehistoric Response 











All site area / model area 0,005 0.0018) ~—~—-0.0022| 0.0029 


nventory area (m*) 218014896.00 54157048.00 115107648.00 387279592.00 
nventory area (km*) ; 218.04 : 
% Inventory area 96.29% 13.98% 100.00% 


o 
oo 


i) anni 2] tow. [Jota 
Model area (m*) —_ 4000722110.04 3066087627.88 §720341377.20| 13787151115.12 
Model area (km*) 4000.72 3066.09 6720.34 13787.15 
% Model area 29.02% 22.24% 48.74% 100.00% 
All sites area (m*) | 19810194.00 5458229.50 15062449.00 40330872.50 
All sites area (km*) 19.81 40.33 
49.12% 13.53% 37.35% 100.00% 
se 2 


‘a 
a 


29,72% 
% Model area inventoried 5.459 1.71% 2.81% 


Oe 






Spring/Steptoe Analytic Unit Model Summary Prehistoric Composite 
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Note; Total area may vary between response and composite analysis due to grid variation within the vegetation evidential theme. 





Figure 5.37 Spring/Steptoe Valley Analytic Unit Composite Probability - Prehistoric 
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Figure 5.38 Spring/Steptoe Valley Analytic Unit - Inventories and Historic Sites 
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Table 5.21 
Spring/Steptoe Valley Analytic Unit Historic Evidential Theme Weights/Chi Square 
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Figure 5.39 Spring/Steptoe Valley Analytic Unit Predictive Pattern - Roads 
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corresponding chi-square value for water is for the buffered area greater than 1000 meters 
from water. (Figure 5.40) Areas within 200 meters of roads and more than 1000 meters 
from water were selected as most predictive for the Spring/Steptoe Valley historic response 
theme. 


Historic Predictive Response 


Three apparent breaks are evident in the normalized posterior response for Spring/Steptoe 
Valley analytic unit historic themes. Breaks lie at 0.006 and 0.0025, with prior probabilities 
set at 0.0026. (Table 5.22) (Figure 5.41) Highest probabilities occur within 200 meters of 
roads, or where proximity to roads and distance to water is greater than 1000 meters 
intersect. No training points fall within a very small zone defined as medium probability. The 
resulting probability map reflects high and low probabilities. (Figure 5.42) The summary 
table (Table 5.23) indicates that 26% of the analytic unit is characterized by high probability 
and 78% of all sites fall within that area. Conversely, the low probability zone covers 75% of 
the area and contains less than 22% of the site area. 


Reclassification of the response theme creates a medium probability zone for the analytic 
unit. It consists of the area defined by the area within 200 meters of roads, or any area more 
than 1000 meters from streams and springs. (Figure 5.43) Forty-three percent of the analytic 
unit falls within the medium probability zone and 66% of all sites lie within that area. (Table 
5.23) The extent of low probability area decreases to 51% of the analytic unit and contains 
16% of the all historic sites, while 6% of the area and 19% of the sites fall within the high 
probability zone. Distribution of inventoried sites is slightly higher in high and low 
probability zones but 80% of the sites still fall within combined high and moderate 
probability areas, 


GREAT SALT LAKE ANALYTIC UNIT 
Analytic Unit Description 


The Great Salt Lake sub-region covers approximately 10.2 million acres (16,079 mi’)/ 
4,164,611 hectares (41,646 km’) within southern Idaho, extreme eastern Nevada, and north 
central Utah. {Figure 3.44) Six hydrographic basins comprise the Great Salt Lake sub-region 
within the study area. The majority of hydrographic units contain lakebed deposits derived 
from the relatively recent Lake Gilbert high stand (10,500 B.P.) and current Great Salt Lake 
shorelines. Slightly more than 10 meters separate the modern and prehistoric shoreline. That 
area comprises 18% of the sub-region. (Figure 5.45) Periodic fluctuations of the Great Salt 
Lake create changing environments along lake shorelines. At elevations between 1290 and 
1310 meters shorelines encroach upon steeper alluvial slopes of surrounding mountain 
ranges, effectively eliminating potential river fed marsh areas (Madsen 1982:208). Six 
hydrographic sub-regions fall within the Great Salt Lake analytic unit. (Table 5.24) 
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Figure 5.40 Spring/Steptoe Valley Analytic Unit Predictive Pattern - Water 
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Table 5.22 
Spring/Steptoe Valiey Analytic Unit Historic Response 
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Figure 5.41 Spring/Steptoe Valley Analytic Unit Historic Response 
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Table 5.23 
Spring/Steptoe Analytic Unit Model Summary Historic Response 




















a eee Be Oe ee ee ee ee 
Model area (m?) 3527228706.14 113945.85| 10259808463.16] 13787151115.15 
Model area (km/) 3527.23 10259.81 43787.15 
% Model area 25.58% 0.00% 74.42% 100.00% 
All sites area (m*) | Spg7202,00) (0.00 1654076.63 7541278.63 
All sites area (cr) eel ee ee 
78.07% 0.00% 21.93% 100.00% 
All site area / model area 0.0017 0.0000 0.0002 0.0005 
Inventory area (m*) 157849264.001 —=—=i0.00}_~—«229430336.00] —_ 387279600.00 
Inventory area (km?) 157.85] —_0.00| 229.43 387.28 
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aga 0.00) 1969211. 
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1337560.36 651651.06[ __-1989211.44 
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Tigh (2) | Medium (1) [Low (0) | Total__ 
Model area (m*) 859979840.00} 5939859456.00] 6978998272.00] 13778837568.00 














Model area (km*) 859.98 5939.86 6979.00 13778.84 
% Model area 6.24% 43.11% 50.65% 
All sites area (m”) 1426929.75] 4950314.00 1163290.25 7540534.00 
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Note: Total area may vary between response and composite analysis due to grid variation within the vegetation evidential theme. 








Figure 5.42 Spring/Steptoe Valley Analytic Unit Observed Probability - Historic 
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Figure 5.43 Spring/Steptoe Valley Analytic Unit Composite Probability - Historic 
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Figure 5.44 Great Salt Lake Anal 
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Figure 5.45 Extent of 
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Table 5.24 
Great Salt Lake Analytic Unit Area 









HYDROGRAPHIC UNIT NAME SUBREGION ACRES HECTARES 

Curlew Valley, Idaho, Utah Great Salt Lake 1,257,783 509,009 
; 6 

Great Salt Lake, Utah Great Salt Lake 1,211,149 490,137 





Skull Valley, Utah Great Salt Lake 518,235 209,723 
Southern Great Salt Lake Desert, Nevada, Utah Great Salt Lake 3,525,768 1,426,83 
eG ae a 
TOTAL 10,290,937 4,164,61 


Northern Great Salt Lake Desert, Nevada, Utah {Great Salt Lake 3,007,645 1,217,15 
Rush-Tooele Valleys, Utah Great Salt Lake 770,357 311,754 


(e) 


N 





Curlew Valley 


The Curlew Valley hydrographic unit lies in the northeastern portion of the Great Salt Lake 
analytic unit. The northern half of the sub-region lies within Idaho; the southern half within 
Utah. Several semi-bolsons comprise this hydrographic unit, all of which slope to the 
southwest and drain into the Great Salt Lake. (Figure 5.46) The hydrographic unit is 
relatively mountainous and is bounded by the Pleasantville Hills, Samaria Mountains, and the 
West Hills to the east. The Promontory Mountains and North Promontory Mountains define 
the southern extent of the Hydrographic unit. Portions of the Raft River Mountains, Black 
Pine Mountains extend into the Curlew Valley hydrographic unit along its western extent; the 
Sublett Range and Deep Creek Mountains extend into the hydrographic unit from the north. 
The Hansel Mountains and North Hansel Mountains extend south through the center of the 
hydrographic unit. Curlew Valley and Hansel Valley are the predominant lowland features of 
this hydrographic unit. Deep Creek, flowing through the upper portion of Curlew Valley is 
the dominant hydrographic feature. The Curlew National Grasslands lie in the upper portion 
of Curlew Valley, where the Sublett Range, Deep Creek Mountains, Pleasantville Hills, and 
North Hansel Mountains merge to form a narrow, well-watered basin. 


Elevations within the Curlew Valley hydrographic unit range from 2429 meters in the Deep 
Creek Mountains to 1285 meters at the Great Salt Lake. Curlew Valley averages 1400 meters 
across its broad southern extent. The upper narrower portion lies at approximately 1580 
meters. Scattered pinyon/juniper woodlands with a sagebrush understory occur on the upper 
slopes of the surrounding and interior mountains. Lowlands range from barren to sparse 
shadscale communities. 


Wetlands are common along the southern periphery of the Curlew hydrographic unit. The 
Bear River National Wetlands extends along the eastern side of the Promontory Mountains 
and the Great Salt Lake. Rozel Flat lies west of the Promontory Mountains and the 
Locomotive Springs State Wildlife Management Area occurs at the delta of Deep Creek and 
the Great Salt Lake. 


Northern Great Salt Lake Desert 


The Northern Great Salt Lake Desert hydrographic unit encompasses the northern half of the 
Great Salt Lake Desert. (Figure 5.47) Only the extreme western edge of the hydrographic 
unit lies within Nevada. Its eastern edge borders the Great Salt Lake while Interstate 80 
arbitrarily bound the southern boundary. The Pilot Range, Goose Creek Mountains and Raft 
River Mountains define the western and northern periphery, respectively. The Leppy Hills lie 
in the southwest corner of the hydrographic unit. Elevations range from 1285 meters on the 
desert floor to 2600 meters in the Pilot Range and 2598 meters at Ingham Peak in the Grouse 
Creek Mountains. Several small ranges lie scattered about the northern Great Salt Lake 
Desert rising as high as 2300 meters. The 1295 meter shoreline of Lake Gilbert (10,500 B.P.) 
roughly defines the edge of the Great Salt Lake Desert sand sheet. The Gilbert shoreline and 
others marking the Lake Bonneville recession are visible along the western ranges and 
mountain “islands” throughout the hydrographic basin. 
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Figure 5. 47 Great Salt Lake Analytic Unit - Northern Great Salt Lake Eucmeachicl Unit 
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Surrounding mountain slopes drain south and east into the Great Salt Lake Desert. Grouse 
Creek Valley and Tecoma Valley provide most consistent drainage systems but both e 
terminate at the edge of the desert. Lowlands along the western edge of the Grassy 

Mountains sustain a viable marsh environment. 


Scattered pinyon/juniper woodlands occur in uplands of the highest interior mountains and 
along the bordering western ranges. As elevation decreases, sagebrush gives way to saltbrush 
communities while most of the bottomlands are barren. 


Great Salt Lake 


The Great Salt Lake hydrographic unit lies wholly within the current extent of the Great Salt 
Lake. (Figure 5.48) Mean elevation for the lake during September 1984 was 1282 meters. 
Land along the periphery of the Great Salt Lake hydrographic unit, if present, consist of 
sandy beach or salt flat. The Bear River and Farmington wetlands border the hydrographic 
unit, but lie outside of its boundaries. Three islands Firemans Island, Antelope Island, and 
Carrington Island, are prominent topographic features in the southern part of the lake. The 
Promontory Mountains form a peninsula within the north-central portion of the hydrographic 
unit and vegetation is sparse to barren. 


Rush-Tooele Valleys 


Rush and Tooele valleys are typical of the north-south trending valleys commonly associated 

with the Great Basin. This hydrographic unit lies south of the Great Salt Lake. (Figure 5.49) 
and consists of the Tooele Valley, a broad open flat sloping northward into the Great Salt ® 
Lake, and Rush Valley, a larger enclosed basin to the south. The hydrographic unit is 

bounded by the Stansbury and Onaou mountains to the west, the Sheeprock Mountains and 

West Tintic Mountains in the south and the Oquirra Mountains to the east. South Mountain 

(2011 meters) divides Tooele and Rush Valleys. Deseret Peak (3362 meters) in the Stansbury 
Mountains and Flat Top Mountain (3237 meters) in the Quirra Mountains provide the highest 

relief along the hydrographic unit boundary, The Tooele Valley continues sloping northward 

from South Mountain with elevations ranging from 1600 meters to 1285 meters at the Great 

Salt Lake. Mud flats and sand sheets dominate the northern portion of the Tooele Valley as it 

juts into the Great Salt Lake. Stansbury Island is a prominent peninsula at the extreme 

northern end of the valley. 


Hydrologically, Rush-Tooele Valley is characterized by steep, well-watered canyons 
draining into the valley floor from the surrounding ranges. Small wetlands and ponds lie at 
the 1520 meter elevation below South Mountain in the northern part of Rush Valley. 
Wetlands also lie at the north end of the Stansbury Mountains and several sloughs grade into 
the mud flats at the north end of Tooele Valley. Vegetation ranges from limber pine at 
highest elevations, pinyon/juniper woodland on slopes above mountain pediments, to barren 
mud flats at lowest elevations. 
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Figure 5.48 Great Salt Lake Analytic Unit - 
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Skull Valley 


The Skull Valley hydrographic unit lies west of Rush-Tooele Valleys, with the crest of the 
Stansbury Range as a common boundary. (Figure 5.50) The Cedar Mountains rising to 2300 
meters, define the hydrographic units western extent, while the Sheeprock Mountains and 
Davis Mountain trend southeasterly to form a southern boundary. The Lakeside Mountains 
form a partial northern boundary. Elevation of the valley floor ranges from 1525 meters in 
the south to 1285 meters in the north where it enters the Great Salt Lake. Extensive mud flats 
dominate the valley floor below 1300 meters in the northern half of Skull Valley. 


Deep canyons along the west slope of the Stansbury Range provide substantia) hydrologic 
inflow, sustaining drainages that eventually flow through the mud flats to the Great Salt 
Lake. Less competent drainages in the Cedar Mountains characterize the hydrologic regime 
along the valleys drier west side. Vegetation ranges from barren on the mud flats to desert 
shrub on the valley floor with pinyon/juniper and sagebrush on the mountain slopes. 


Southern Great Salt Lake Desert 


This hydrographic unit is the southem extension of the Northern Great Salt Lake Desert 
hydrographic unit. (Figure 5.51) It shares similar characteristics; dry mud flats and a sand 
sheet comprise most of the unit, but fewer “islands” occur within its interior. The Goshute 
Mountains, Ferber Hills, and Deep Creek Range form the western boundary of the Southern 
Great Salt Lake Desert. The Leppy Hills and Danger Cave mark the extreme northwest 
corner of the hydrographic unit. To the east, the hydrographic unit boundary is shared with 
the ranges bordering Skull Valley. Several low, north-trending ranges extend into the Great 
Salt Lake Desert, creating the hydrographic unit’s southern boundary. Deep Creek Valley, 
Snake Valley, Fish Springs Flat, and Dugway Valley lie between these southern ranges and 
drain northward into the desert. White Horse Flat and related badlands lie between the 
Goshute Mountains and Ferber Hills. Highest elevations occur within the Goshute 
Mountains, with Goshute Peak rising to 2929 meters. Southern valleys slope northward with 
highest elevations between 1600 and 1550 meters, The desert floor where at the boundary 
with Northern Great Salt Lake Desert is 1285 meters. Wildcat Mountain and Granite Peak 
(2154 meters) are “island” features within the hydrographic unit. 


Intermittent streams originating in the surrounding mountains provide water flow into the 
hydrographic basin. Sustainable wetlands occur at the north end of White Horse Flat where 
Felt Wash, originating in the Goshute Mountains, enters the Great Salt Lake Desert mud 
flats. That marsh lies at 1290 meters. Numerous springs feed the lowlands of Fish Springs 
Flat at an elevation of 1309 meters along the eastern terminus of the Fish Springs Range and 
the Fish Springs Wash delta. 


Vegetation is typical of the Great Salt Lake Desert. Lowest elevations are barren mud flats 


and sand sheets, grading to desert shrub communities as elevation rises from the desert floor. 
Pinyon/juniper uplands grade sagebrush communities along alluvial fans. 
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Figure 5.50 Great Salt Lake Analytic Unit - Skull Valley Hydrographic Unit 
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Figure 5.51 Great Salt Lake Analytic Unit - Southem Great Salt Lake Hydrographic Unit 
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Analytic Results 
Prehistoric Evidential Themes 


Approximately 1062 square kilometers, 2,5% of the total area, were inventoried within the 
Great Salt Lake analytic unit. (Table 5.25) (Figure 5.52) Three hundred eleven sites from 
those inventories were considered in the weighted analysis. One thousand one hundred 
sixteen sites are reported for the entire analytic unit. 


Within all analytic classes, only Great Basin pine, Juniper steppe and Wet grassland classes 
within the vegetation evidential theme have not been inventoried. (Table 5.26) The total 
extent of the missing class area is 202 square kilometers (0.4% of the total analytic area). 
Inventoried space in several zones, including chaparral and water in the vegetation theme, 
mountain areas in landform and slopes greater than 15 degrees are under-represented within 
the inventoried sample. When evaluating weighted contrasts, sampling inconsistencies and 
site densities relative to the model area data set as well as inventoried data sets were 
considered. 


Weights of evidence tables identify classes within each evidential theme that lie “inside” of 
the predictive pattern. (Table 5.26) Normalized contrast for juniper/pinyon vegetation class is 
highest when evaluated with all categories of prehistoric sites. Chi-square is also 
significantly high for the class. Barren areas, or those with sparse vegetation, have a 
correspondingly high negative contrast value, indicating a lower than expected probability 
for sites. By contrast, barren areas have the highest contrast for non-sites. (Figure 5.53) 


Contrasts for distance to springs and streams is variable across each different analytic run. 
When all sites are considered, the 200 meter buffer has the highest contrast. When 
inventoried sites are weeded, the 400 meter buffer distance has the highest contrast, while the 
sites within 250 meter cells have highest contrasts in areas greater then 2000 meters. The 
large expanse of desert within the Great Salt Lake analytic unit and peculiarities of the 
weeding process appear to be driving the contrast results. Since the weighted results are 
inconsistent across all analytic runs, distance to water was not included as a predictive theme. 
(Figure 5.54) 


Proximity of sites to wetlands, on the other hand, uniformly identifies the 0 to 1000 meter 
buffer as a reliably predictive class. (Figure 5.55} In the three analytic runs with sites, the 
1000 meter buffer exhibits the highest contrast values. Corresponding negative values are 
present in the non-site analysis. Site location more than 5000 meters from potential wetlands 
is inconsistently identified in the weights tables. Relatively high contrasts in weeded all site 
and inventoried site analysis, likely reflect an upland adaptation within the analytic unit. 


Like distance to streams and springs, analytic runs for slope are less than conclusive. 
Contrast values for slopes 15 to 30 degrees are based upon a relatively high frequency of 
sites within a slope class that accounts for less than 1% of the analytic unit. Those sites are 
most likely rockshelters. Contrast for slope between 5 and 15 degrees is also relatively high, 
but chi-square calculations suggest that the distribution of sites within that class is normal. 
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Great Sait Lake Analytic Unit inventory Summary 
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Figure 5.53 Great Salt Lake Analytic Unit Predictive Pattern - Vegetation 
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Figure 5.54 Great Salt Lake Analytic Unit Predictive Pattern - Streams and Springs 
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Figure 5.55 Great Salt Lake Analytic Unit Predictive Pattern - Potential Wetland 
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Most of the sites within the analytic unit lie on slopes between 0 and 5 degrees, but 
proportionally, that frequency lies outside of the expected distribution. While one would be 
inclined to accept the 15 to 30 degree slopes as predictive for the theme, non-site analysis 
suggests that it is likewise predictive for non-sites. As a result, slope was not selected as a 
predictive theme in the Salt Lake analytic unit. 


Landform as an evidential theme provides stronger relational contrasts than incremental 
slope. While most of the sites occur within flats along the valley floors, the piedmont is 
consistently characterized by higher than expected site frequencies. Chi-square also confirms 
a non-random distribution. The piedmont as a landform class subsumes a portion of the 
slopes within a 5 to 30 degrees range and was selected as a predictive evidential class. 
(Figure 5.56) 


Prehistoric Predictive Response 


After identifying evidential theme classes that lie inside the predictive pattern, a response 
theme was calculated in order to compile a probability map based upon the likelihood of 
encountering a site within the aggregated evidential themes. (Figure 5.57) Normalized 
posterior probabilities were used as a means to evaluate tabular results. Observed breaks are 
apparent at posterior probabilities of 0.016 and 0.004. (Table 5.27) (Figure 5.58) Highest 
probabilities for encountering sites occur in areas within 1000 meters of potential wetlands 
and on piedmont slopes, or a within a combination of proximity to potential wetland, 
piedmont slopes, and juniper/pinyon vegetation zones. Probabilities decrease in areas 
characterized by the presence of a predictive single evidential theme. When no predictive 
classes are present, posterior probabilities fall below theO .005 critical prior probability value, 
and predictive probabilities are lowest. 


Results of the probability model were analyzed in a spatial context in order to validate model 
results. (Table 5.28) Extent of the sensitivity areas, those within both the entire analytic unit 
and inventoried portions of the analytic unit, were contrasted with actual areal extent of the 
sites within each sensitivity zone. Highest ratios of site area to sensitivity area should fall 
within zones of highest probability if the model is accurate. Summary tables show that the 
areal density of all sites and inventoried sites are indeed highest within areas of high to 
medium sensitivity. While areas of low sensitivity comprise two-thirds of the analytic unit 
they consistently maintain the lowest values of site to total area. (Figure 5.59) 


The response theme calculated with Spatial Data Modeler accurately grouped the 
intersection of predictive classes into probability zones. As a result, recalculating the 
response using Spatial Analyst® produce similar results. Variation in probability and site 
areas change less than 0.5%. 


Historic Evidential Themes 
Two hundred three historic sites are reported within the Great Salt Lake analytic unit. Within 


the 1062 square kilometers subset of inventories greater than 640 acres in extent, 61 sites are 
considered for analysis. (Table 5.25) (Figure 5.60) Distance to existing roads and water 


V-15 








Figure 5.56 Great Salt Lake Analytic Unit Predictive Pattern - Landform 
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Figure 5.57 Great Salt Lake Analytic Unit Observed Probability - Prehistoric 
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Table 5.27 
Great Salt Lake Analytic Unit Response 
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Figure 5.58 Great Salt Lake Analytic Unit Observed Breaks 
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Table 5.28 
Great Salt Lake Analytic Unit Model Summary Prehistoric Response 


Sr a ee ls a a 
Model area (m*) 1588656090.13] 12304292627.13] 27753106351.78] 4164605506 
Mode! area (km”*) 1588.66 12304.29 27753,.41 41646.0 
% Model area 66.64% 100.00% 
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All sites area (km*) 11.97 
11,78% 43,10% 45.12% 400.00%| 
All Site area/modelarea | 0.0021 0.0010 0.0007] 
47948860.00] —-363288416.00| —_ 650945792.00/ _ 1062183068.0 

nventory area (km’) 363,29 650.95 __ 1062.18 
4.51% 34.20% 61.28% 100.00% 
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_ Model area inventoried 3.02% 2.35% 2.55% 


000 
% 
nventory sites area (m*) | __621116.56| 3849582,25 2648310.00 7119008.81 
mveniorysites area’) | __oea|——awal =a] a2 
% Inventory site area 8.72% 54.07% 37.20% 100.00% 
nv site area / inv area 0.0130 0.0106 0.0041 0.0067 















| C—~—~—SSSCCsHiigh (3-2) | Medium (1) | Low(0) | 
Model area (m*) 1599121920.00| 12245988352.00] 27664097280.00] 41509207552.00 

























Model area (km’) 1599.12 42245.99 2766410 41509.21 
% Model area 3.85% 29.50% 66.65% 100,00%| 
All sites area (m’) 3270916.00]  11951653.00]  12302428,00|] _—- 27524997.00 
All sites area (km*) Paar Eee ee 12.30 27.52 
43,42% 44.70% 100.00% 
All site area / model area 0.0010 0.0004 0.0007 
inventory area (m’) 47983864.00}  362866880.00{  650192896.00| _1061043640.00 
Inventory area (km?) 47.98 362.87} «650.19 1061.04 
4.52% 34.20% 61.28% 100.00% 
% Model area inventoried 3.00% 2.96% 2.35% 
inventory sites area (m') 621116.56 3849582.25 2648310.00 7119008.81 
inventorysites area’) | _ oe —ssasl SCs 
% Inventory site area 8.72% 54.07% 37.20% 00.00% 
Inv site area / inv area 0.0106 0.0041 0.0067 


Note: Total area may vary between response and composite analysis due fo grid variation within the 
vegetation evidential thame. 
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Figure 5.59 Great Salt Lake Analytic Unit Composite Probability - Prehistoric 
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sources were considered as predictive evidential themes for historic resources. (Table 5.29) 
Roads were buffered at 200 meter intervals to 1000 meters, and water sources at 200, 400 t 
and 1000 meter intervals. Like prehistoric sites, weights were calculated for using all sites, 

inventoried sites, and inventoried site areas within a 250 meter inventoried grid. Within the 

250 meter grid, non-site weights were also calculated for comparison with weighted site 

results and calculation of chi-square. 


Within the road evidential theme, a buffered distance of 200 meters consistently revealed the 
highest positive contrast, with a corresponding negative contrast for non-sites. (Figure 5.61} 
Chi-square results for this class fail to meet the desired threshold for non-random 
distribution. Highest contrasts are reflected within the 600-800 meter buffer for site areas 
within the 250 meter grid, but a high positive contrast in non-site weights within the same 
class suggests that the occurrence of non-sites is highly probable for that buffered area. With 
the conflicting results, no class within roads was clearly predictive for a positive contrast. 


Within inventoried areas, contrasts for distance to water showed increasingly high probability 
as distance from the water source increased from 0 to 1000 meters. Chi-square results show a 
correspondingly positive relationship. Negative contrasts for non-site areas corroborate the 
probability of encountering sites within an ascending radius of water sources. (Figure 5.62) 


Historic Predictive Response 


Response themes cannot be built without a minimum of at least two evidential themes. The 

water buffer alone could have been used as a predictive mask, but further analysis of the 

weights tables showed a good correlation for sites not occurring more than 1000 meters from YY | 
roads, or more than 1000 meters from water sources. The highest negative contrasts for each 

theme was selected as “inside” the pattern and used to calculate a response theme. Since the 

resulting response table calculates weights and probabilities for combined classes, the only 

difference between selecting positive or negative contrasts is that the intersection of 

evidential theme classes appears at the bottom of the table, associated with the lowest 

posterior probabilities. 


Three observed breaks are evident within posterior probabilities generated for historic 
evidential themes in the Great Salt Lake analytic unit. (Table 5.30) (Figure 5.63) Breaks 
occur at 0.002 and 0.0007, with prior probabilities set at 0.0003. No areas within high 
probability fall within buffers lying further than 1000 meters from water or roads. Moderate 
probability areas lie within 1000 meters of roads, but occasionally more than 1000 meters 
from water, and low probability zones always occur more than 1000 meters distant from 
roads, or more than 1000 meters from roads and water. 


The corresponding sensitivity map (Figure 5.64) and summary table (Table 5.31) shows that 
areas of highest and medium probability include well over 80% of all sites and inventoried 
sites by area. Slightly less than 50% of the analytic unit falls within the low probability area, 
yet less than 10% of the historic site areas occur within that zone. 











Great Salt Lake Analytic Unit Historic Evidential Theme Weights/Chi-Square 
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Figure 5.61 Great Salt Lake Analytic Unit Predictive Pattern - Roads 
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Figure 5.62 Great Salt Lake Analytic Unit Predictive Pattern - Water 
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Table 5.30 
Great Sait Lake Analytic Unit Historic Response 
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Figure 5.63 Great Salt Lake Analytic Unit Historic Response 
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Figures 5.64 Great Salt Lake Analytic Unit Observed Probability - Historic 
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Table 5.31 
Great Salt Lake Analytic Unit Model Summary Historic Response 


























Model area (m*) 15518336331.13| 5401380212.77| 20726338525.13| 41646055069.03 
Model area (km*) 15518,34 5401.38 20726.34 
37.26% 49.77% 10.00% 
Alt sites area (m*) 8448377.00 2772682.25 948804.00| _ 12169863.25 
All sites area (km*} ee ae eee) 12.17 
69.42% 22.78% 7.80% 400.00% 
0.0005 0.0005 0.0000 90,0003 
Inventory Area (m*} 377139904.00]  212271424.00] _472771744.00| _1062183072.00 
Inventory Area (km*) 212,27 1062.18 
35.51% 19.98% 10.00% 
% Model area inventoried 2.43% 3.93% 2.28% 
‘Inventory sites area (m*) §553557.00 —-452059.66 266618.38 6272235,.03 
(nventory sites area(km’) [| 5.55 0.45 
% Inventory site area 98.54% 7.21% 4.25% 00.00% 
Inv site area / inv area a aa es 0.0006 ___0.0059 





Great Salt Lake Analytic Unit Model Summary Historic Composite 
































i a) High @).. | Medium (4) | Low (0). | Total. 
Model area (m*) 415518270464.00| 14305559552.00] 11821988864.00| 41645818880.00 
Model area (km’) 15518.27 14305.56 11821.99 41645.82 
37.26% 34.35% 28.39% 4100.00% 
All sites area (m’) 8448377.00 3560621.50 160864.72 12169863.22 
Al sites area (km?) ames ee eel eel 12.17 
69.42% 29.26% 4.32% 100.00% 
All site area / model area 0.0005 0.0002 0.0000 0.0003 
Inventory area (m”) _377139904.00] _344971424.00] _ 340071008.00] _1062182336.00 
Inventory area (km*) 377.14 344.97 340.07 1062.18 
% Inventory area 35.51% 32.48%] «32.02% 400.00% 
% Model area inventoried 2.43% 2.41% 
inventory sites area (m’) 5653557.00 685164.56 33513.48 6272235.05 
Inventory sites area (km*) Pee ee | ee a 
% Inventory site area 88.54% 10.92% 100.00% 
9.0020 0.0001 0.0059 








Recalculating the historic probability by intersection of predictive classes redefines medium 
and low probability zones. The medium probability area is expanded to 34% of the analytic 
unit with 29% of all sites present, while the low probability zone decreased to 28% of the 
analytic unit and includes slightly more than 1% of all sites. Within inventoried areas, less 
than 1% of the sites fall within the low probability zone. (Figure 5.65) 


UPPER SNAKE ANALYTIC UNIT 
Analytic Unit Description 


The Upper Snake sub-region covers approximately 3.0 million acres (4801 mi’)/1.2 million 
hectares (12,435 km”) within southern Idaho, northeastern Nevada, and northwestern Utah 
(Figure 5.66) or 15% of the GBRI study area. Three hydrographic units, Salmon Falls, 
Goose, and Raft comprise the analytical portion of the sub-region. (Table 5.32) All three 
drain in a northeasterly direction towards the Snake River. Complex, dendritic drainage 
patterns dominate the Upper Snake sub-region. 


Salmon Falls 


The Salmon Falls hydrographic unit lies in the westernmost portion of the Upper Snake 
analytic unit. Salmon Falls Creek and its tributaries is the dominant hydrologic feature of the 
hydrographic unit. (Figure 5.67} Elevations range from 2631 meters at Ellen D Mountain, 
and 2410 meters at Middle Stack Mountain near Contact, Nevada, in the southern portion of 
the hydrographic unit to 900 meters at the confluence of Salmon Falls Creek and the Snake 
River in the northern portion of the hydrographic unit. Major physiographic feature include 
the O’Neil and Shoshone Basins, Antelope Pocket and Browns Bench, a major obsidian 
source, all within the southern half of the unit. Vegetation is primarily sagebrush with some 
pinyon/juniper woodland. Topography becomes more subdued progressing northward 
through the Hydrographic unit. Higher mountains give way to low ridges and dissected basalt 
plateaus. 


Goose 


The Goose hydrographic unit lies in the central portion of the Upper Snake analytic unit, 
covering portions of Idaho, Nevada, and Utah. (Figure 5.68} Goose Creek and its tributaries 
dominate the hydrology of this hydrographic unit. Like Salmon Falls Creek, it drains 
northward towards the Snake River. Lowest elevations (1290 meters) occur in agricultural 
lands near the Snake River. Monument Peak (2454 meters) lies in the uplands within the 
mountainous, west central portion of the hydrographic unit. The Sawtooth National Forest 
administers most of this area. To the south, low hills and ridges characterize the hydrographic 
unit, while the northern one-third is relatively flat agricultural lands. Deadman Ridge and 
Middle Mountain flank respective western and eastern edges of the unit, while Big Draw and 
Cedar Mountain Draw lie in the south. Sagebrush dominates the landscape outside of 
agricultural areas, pinyon/juniper woodlands are found in the steeper uplands. 
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Figure 5.65 Great Salt Lake Analytic Unit Composite Probability - Historic 
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Table 5.32 
Upper Snake Analytic Unit Area 
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Figure 5.68 Upper Snake Analytic Unit - Goose Hydrographic Unit 
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Raft 


Along the eastern side of the Upper Snake analytical unit, the Raft River and its tributaries 
create a major hydrologic feature. (Figure 5.69) The hydrographic unit lies within Idaho and 
Utah. Several ranges including the Jim Sage Mountains, Alison Mountains, Black Pine 
Mountains, and Middle Mountain bound this horseshoe-shaped basin. The Raft River 
Mountains provide a topographic divide between the Upper Snake sub-region and the Great 
Salt Lake sub-region to the south. Highest elevations occur at Cache Peak (3151 meters) and 
Mount Independence (3033 meters) along the western edge of the hydrographic unit. The 
Upper Raft River Valley, Junction Valley and the Holt Basin lie in the southern portion of 
the hydrographic unit. Basins in the south average approximately 1700 meters, while 
agricultural lands in the Raft River Valley lie at 1285 meters near the confluence of the Raft 
and Snake Rivers. Juniper and pinyon dominate higher elevations across the hydrographic 
unit. Sagebrush is the dominant non-cultivated plant community. 


Analytic Results 
Prehistoric Evidential Themes 


Since a major portion of the Upper Snake analytic unit contained no spatial inventory data, 
statistical relationships between site and non-site components could not be evaluated. Site 
location data was derived primarily from Bureau of Land Management data sets. Forest 
Service and other agency lands within the analytic unit were not included in the analysis. 


One thousand six hundred seventy-five prehistoric sites have been recorded on 8085 square ( 
kilometers of BLM land within the Upper Snake analytic unit. (Table 5.33) Weights tables 

for evidential themes were compiled using all weeded sites so that only one training point 

would occur within each analytic cell. Weeding reduced the total number of training points 

by approximately 20%. In each of the evidential themes the class with the highest positive 

contrast was selected as most predictive, with the remainder falling outside of the pattern 

regardless of whether contrast were positive or negative. (Table 5.34) 


Within the vegetation evidential theme, the juniper steppe class and sagebrush zone have a 
positive contrast, while juniper/pinyon is least predictive for sites. Juniper steppe was 
considered inside the pattern since its contrast and positive weights were highest. The area of 
juniper is relatively small, 5.6% of the total analytic unit and weeded sites account for 10.6% 
of the total. (Figure 5.70) 


The only positive contrast for distance from streams and springs is within the 0 to 200 meter 
buffer. Areas more than 1000 meters from water courses exhibit the highest negative 
contrast. Over 56% of the weeded sites are located within the 0 to 200 meter buffer. (Figure 
5.71) 


Areas within 1000 meters of potential wetlands, while relatively small (10.3% of the analytic 
unit) also have the highest contrast. Areas lying more than 5000 meters from potential 
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Upper Snake Analytic Unit Site Summary 
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Figure 5.69 Upper Snake Analytic Unit : Raft Hydrographic Unit 
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Figure 5.70 Upper Snake Analytic Unit Predictive Pattern - Vegetation 
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Figure 5.71 Upper Snake Analytic Unit Predictive Pattern - Streams and Springs 
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wetlands cover the greatest proportion of the analytic unit. They are also moderately 
predictive, but considered outside of the probability pattern. (Figure 5.72) 


Slope and landform exhibit contrasting predictive results. Two percent of the analytic unit 
lies on slopes between 15 and 30 degrees, but almost 6% of the weeded sites occur within 
that slope class. Most of the area and most of the sites fall within the 0 to 5 degree slope 
class. It has a high negative contrast. (Figure 5.73) Within landform, however, areas along 
the basin or valley floors have the highest contrast. Landform classes are evenly distributed 
across the analytic unit, but piedmont and mountain both have negative predictive contrasts. 
(Figure 5.74) 


o 


Prehistoric Predictive Response 


Snake response themes were mun using most predictive classes for each evidential theme. 
Three distinct breaks occur within the normalized posterior probability values. Highest 
probabilities range between 0.24 and 0.17, medium probability ranges between 0.17 and 
0.058, while lowest probabilities fall between 0.058 and 0,021. Prior probability was set at 
0.045. (Table5.35) (Figure 5.75) 


The resulting table shows that less than 1% of the analytic unit falls within the high 
probability zone and a similar percentage of sites are associated with that area while more 
than one-half of the sites fall within the low probability zone which extends over 66% of the 
analytic unit. (Table 5.36) Site densities within relatively large cells within the low 
probability zone (values 8 and 7) (Table 5.35) likely cause the skewed response, (Figure 
5.76) 


By calculating an intersection of predictive themes, a better correlation between site density 
and probability is gained. The total area of high probability is expanded to almost 19% of the 
analytic unit and it contains nearly 24% of the sites, while the low probability zone is 
decreased to 12% of the analytic unit and includes 8% of the sites. Most of the analytic unit 
lies within areas of medium probability (70%) and most of the site areal falls within that 
zone. (Table 5.36) (Figure 5.77) 


Historic Evidential Themes 


One hundred nineteen historic sites were recorded within the 8085 square kilometers of land 
managed by BLM. (Table 5.32) Within both the road and water evidential themes considered 
for historic resources, the buffered class between 0 and 200 meters is most predictive. 
Contrasts within roads decline consistently with each increasing buffer. (Figure 5.78) A 
similar, but not as striking decline occurs with buffered distance to water courses. (Figure 
5.79) (Table 5.37) 


Historic Predictive Response 


Three breaks are evident in the historic response theme when the area within 200 meters of 
roads and water are selected as predictive classes. Breaks occur at normalized posterior 
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Figure 5.72 Upper Snake Analytic Unit Predictive Pattern - Potential Wetland 
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Figure 5.73 Upper Snake Analytic Unit Predictive Pattern - Slope 
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Figure 5.74 Upper Snake Analytic Unit Predictive Pattern - Landform 
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Figure 5.75 Upper Snake Analytic Unit Observed Breaks 
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Table 5.36 
Snake Analytic Unit Model Summary Prehistoric Response 

















Model area (m*) _ 66845996.19| 2688005624.89] 6330581114.90| 8085432735.98 
Model area (km?) 66.85 2688.01 5330.58 8085.43 
0.83%| = 33.25%] ———=—«65.93%| 10.00% 
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All sites area (km!) __17.86 
16.10%) sa.00%| 100.00 
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Snake Analytic Unit Model Summary Prehistoric Composite 
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Figure 5.76 Upper Snake Analytic Unit Observed Probability - Prehistoric 
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Figure 5.77 Upper Snake Analytic Unit Composite Probability - Prehistoric 
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Figure 5.78 Upper Snake Analytic Unit Predictive Pattern - Roads 
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Figure 5.79 Upper Snake Analytic Unit Predictive Pattern - Water 
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Upper Snake Analytic Unit Historic Evidential Theme Weights 
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probabilities of 0.006 and 0.0027, the prior probability was set at 0.0033. (Table 5.38) 
(Figure 5.80) Summary tables show that the high probability area has the smallest extent 
(less than 10% of the analytic unit) while 77% of the analytic unit falls within the low 
probability zone and 13% falls within medium probability. (Figure 5.81) Sites are evenly 
distributed across probability zones with approximately one-third of the sites within each 
zone. (Table 5.39) Proximity to water but not roads, and distances greater than 200 meters 
for either evidential theme, fall within the low probability zone. Again, the number of 
training points relative to low probability area biases the response pattern. 


Summary calculation of the predictive classes creates a modified response with better site to 
area ratios between probability zones. Training points associated with distance to water but 
not roads are included in the medium probability zone with the summary calculation, The 
composite summary produces a 25% reduction to the areas of the low probability zone and 
corresponding increase in the extent of the medium zone. (Table 5.39) One-third of the sites 
still remain in the high probability zone, while almost one-half fall within zones of moderate 
probability. The remainder of sites lie within the low probability zone. (Figure 5.82) 
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Figure 5.80 Upper Snake Analytic Unit Observed Breaks- Historic 
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Table 5.38 
Upper Snake Analytic Unit Historic Response 
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Figure 5.81 Upper Snake Analytic Unit Observed Probability - Historic 





[_] Snake Analytic Unit 
2 Historic Sites 


Probability 10 0 10 20 Miles 
Low 


ees Medium 20 0 20 40 Kilometers 
ee High a 
[_] No Data 


Table 5.39 
Upper Snake Analytic Unit Model Summary Historic Response 
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Figure 5.82 Upper Snake Analytic Unit Composite Probability - Historic 
["_] Snake Analytic Unit 
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VI. DISCUSSION 


The GBRI cultural resources planning model provides a statistically useful indicator for 
predicting the likelihood of cultural resources on a landscape level. The statistical 
technique of Weights of Evidence provides an adequate means to evaluate the 
distribution of sites within the chosen evidential themes, but some caution is required due 
to biases resulting from site to unit area ratios. That bias hold through calculations of 
predictive responses. Intersecting themes provided the most reliable means of identifying 
probability. The probability of encountering cultural resources is highest in areas where 
multiple predictive evidential themes intersect, while the likelihood of encountering 
cultural resources lessens as fewer predictive themes are encountered. 







Table 6.1. Summary of weighting factors for prehistoric site 
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*High slope values are always a negative factor — 


Table 6.1 identifies predictive classes for prehistoric cultural resources within each 
evidential theme for each analytic unit within the study area. Each predictive class 
identifies a landscape element as a potentially predictive surface, and combinations of 
predictive surfaces increase the probability of encountering a cultural resource within that 
area. For example, within the Pilot/Thousand Springs Analytic Unit, activities within the 
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Table 6.2 
Distribution of Prehistoric Site Area by Probability Zone 























Ce a OC ea eae 1): eas 
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Note: Total area may vary between prehistoric and historic composite analysis due fo grid vanation within the 
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Distribution of Historic Site Area by Probability Zone 
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Table 6.3 Spring/Steptoe Model Summary and Additional Data 


Model Summary 
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desert shrub vegetation zone, lying upon piedmont landforms, and within 1000 meters of 
potential wetland areas would have the highest probability of encountering sites, while 
sites less likely to occur in areas where those conditions are not met. With the possible 
exception of potential wetlands, the evidential classes are readily identifiable at the field 
level and provide a basis for evaluating probability of encountering cultural resources. 


Historic resources are much easier to predict — they are nearly always within 1000 meters 
of perennial water sources and (not surprisingly) within 500 meters of roads. 


An analysis of composite probabilities for the study area as a whole quantifies the 
effectiveness of the model. For prehistoric sites, high probability zones cover 17% of the 
entire study and contain 33% of the all site areas, while low probability zones extend over 
60% of the study area, but contain only 25% of the sites. (Table 6.2; Figure 6.1) The 
remaining 41% of the site area falls within 46% of the study area. Historic sites and 
probability zones exhibit a similar pattern. (Table 6.2; Figure 6.2) Twenty-five percent 
of the study area contains 50% of the historic site area while 40% of the study area and 
8% of the historic site area falls within the low probability zone. 


From its inception, the GBRI cultural resources probability model functions was 
conceived as a pattern recognition tool rather than as an explanatory model relating to 
human adaptive response. Buffers within evidential themes were chosen since they 
represent potential foraging radii, or in the case of slope, habitable ground, but the results 
are never synthesized to suggest a causal relationship. Evidential themes provide only a 
recognizable landscape layer that can be contrasted against site density patterns. 


Since the model is based upon pattern recognition, subsequent inventories and new site 
data may provide subtle, or in some cases dramatic changes to the distributional patterns. 
Certain classes within an evidential theme may have been inadequately sampled during 
previous investigations, or sites poorly reported. Newly acquired data may effectively 
increase both inventoried strata and drive results towards more or less predictable 
distributions. Recently acquired data supports that proposition. 


The of the Spring/Steptoe Valley analytic unit revealed that in addition to several other 
themes, flats, within landform, were moderately predictive. A recent field investigation, 
not included as data in generating the current model, identified 163 sites within 
approximately 20,550 acres (83.16 km’) along the eastern slope of the Egan Range south 
and west of Ely, Nevada. Approximately 4300 acres (17 km’) were previously 
inventoried, resulting in a net increase of 16,250 acres (65.77 km’) of inventory. (Table 
6.3) Most of the inventory was conducted within the piedmont landform, effectively 
increasing the investigated area within that evidential class from 2.7% to 5.7% of the 
analytic unit Additional sites reported within the piedmont, increase site distributions 
within that zone from 20% of all sites (using the model data) to 27% of all sites within 
the analytic unit. Slightly more than 20% of the analytic unit comprises the piedmont 
zone, 
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Figure 6.1 Distribution of Prehistoric Site Area by Probability Zone 
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Figure 6.2 Distribution of Historic Site Area by Probability Zone 
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With the original data, weights tables show that flats are strongly predictive while other 
landforms reveal a negative contrast. (‘Table 6.4) When the new sites are included in the 
analysis, 27% of the sites fall within 20% of the analytic unit characterized as piedmont. 
Resulting contrasts within flats and piedmont are similar — indicating that the piedmont is 
now revealed to be “more sensitive” than before. Flats remain slightly more predictive 
than piedmont. 


The important lesson of the example above is that we can anticipate the model to change 
as more information becomes available. A second survey again helps to illustrate how 
this can occur, especially when initial inventories are small. In Spring/Steptoe, the model 
data contains 65 sites within 77 square kilometers of inventoried piedmont (16% of all of 
the sites in the analytical unit, 20% of the inventoried area within the analytical unit as a 
whole). However, a single new inventory that covered 35 square kilometers of piedmont, 
revealed 102 new sites. With the inclusion of this new data, 29% of all sites fall within 
the piedmont zone, and the piedmont zone comprises 25% of all inventoried ground. The 
statistical significance of the ratio differences (percent of sites: percent of inventory) may 
be questionable, but it illustrates how new inventory will change our picture of specific 
analytical units. 


Regardless of the effect of the new data, the model derived for Spring and Steptoe 
Valleys did adequately predict probability of encountering resources within the new 
inventory. While the spatial extent of sites that were identified during the new 
reconnaissance was not available, each site was buffered to a 2 acre extent so that site 
area per probability zone could be calculated, (Table 6.5) Over 50% of the inventory 
area falls within the medium probability zone and less than 12% falls within areas of high 
probability. Utilizing the derived site extent, 76% of the site areas lie within high to 
moderate probability zones with the remainder falling within the low probability area that 
accounts for 38% of the inventoried extent. 


In addition to providing probability layers useful for long range planning, the model also 
brings together site and inventory information useful for short and long term planning. 
Summary tables and related shapefiles identify the percentage of inventory within each 
analytic umit, and assess the relative densities of site area to cumulative inventory blocks 
and within each analytic unit. They also identify proportional survey coverage within 
specific environmental settings, allowing the cultural resources manager to better assess 
the range of coverage within a resource area. 


Field experience and expert knowledge of the regions within the project area provide the 
best means to verify model results. If regional expertise has intuitively predicted that 
most sites are found within 200 meters of water sources, and the evidential theme reflects 
a similar pattern, then that theme is most likely valid. Likewise, the model may direct 
confirmatory evaluation. If a composite theme is identified as predictive in the model, but 
has never been explored or evaluated by regional experts, subsequent projects can be 
tailored to validate the model's findings. 
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Table 6.5 


Additional Reconnaissance Model Summary (Composite) 
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The IMACS assemblage data provides a useful tool for deciphering cultural patterns and 
compilation of overview information. Unfortunately, the quality of data and its 
completeness are variable, often dependent upon age of the record. The IMACS encoding 
form itself also lacks a level of information that could answer more specific research 
questions. Lithic assembalge characteristics, such as frequency by material type, are not 
preserved in the encoding format and the assessments of lithic stages are inconsistent. 
Older site records pose additional constraints to completeness of the assemblage database 
since IMACS classifications must be derived from narrative descriptions. 


Several problems were also encountered in the creation of a comprehensive assemblage 
database that was compiled from electronic data maintained by three separate entities. 
Administrative data is consistent across the three database used in the analysis, but 
assemblage data varies from complete IMACS encoding to descriptive summaries of the 
cultural assemblage. Consistency in reporting National Register eligibility also varies 
between agency and archive. In some cases, current status is maintained in IMACS 
format within the site database, in others a separate database contains that information. 


Shortcomings of the databases can be overcome in future projects by scaling database 
contents to fit the project goals. A broad based predictive model can be constructed from 
existing assemblage data with minimal effort if research questions are limited in scope. 
Where do we expect prehistoric sites? Where do we expect historic sites? Where do 
National Register sites occur? More detailed synthesis requires mining data from any 
combination of existing electronic data and paper records for completeness and missing 
information. What types of materials comprise the lithic debitage? Are lithic tools 
manufactured from materials available locally? 


Scale of the analytic area should also be adjusted to fit the research questions. Another 
problem encountered with the anthropological analysis was the validity of generalizing 
results to fit such a broad research area. Variations in survey quality, site reporting, and 
archival data over an area in excess of 78,000 square kilometers can only elucidate very 
general patterns. As research questions become more pointed, the research area needs to 
be scaled down, evidential themes refined to be more specific, and site information 
scrutinized to assure validity of the observations. Patterns unique to the Upper Snake 
hydrologic unit may not be valid for the Southern Great Salt Lake Desert. 


Planning Models As Cultural Resource Forecasts 


The GBRI cultural resources model study is phrased as two map layers (distributed as 
Environmental Systems Research Institute (ESRI) GRID format files (Ancillary CD 1). 
One map layer is for prehistoric resources, the other for historic resources. The user of 
these map layers must be well aware of several important aspects of them. 


The map layers are summaries of models, not “known” data. Just as fire managers do not 
really know the accuracy of their fuel regime models until fire actually consumes a spot, 
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we cannot know how accurate the models are currently. The examples above suggest that 
survey bias, differences in reporting styles (sites vs. isolates), and imprecision in the 
baseline data will all contribute to inaccuracy. Below, we discuss long-term strategies for 
coping with these problems; here, we wish to call attention to the nature of the models 
and maps. 


Models and maps are planning tools, not compliance tools. One cannot use the GIS data 
to say an area will be devoid of cultural resources just because it has a LOW value 
associated with it. For these reasons, we prefer to call the GBRI model a planning model 
rather than a predictive model. The maps (paper or electronic), which summarize are 
current planning-level knowledge, are thus forecasts. The simile to meteorology is not 
accidental, for we do not fully understand the system that generated the cultural resources 
we are attempting to forecast. Yet, just as a forecaster can state that a particular weather 
patter is highly likely to yield snow in the Sierra — without necessarily understanding 
why the pattern occurs — the GBRI model can forecast areas of highest and lowest 
likelihood of cultural resources. If one thinks of the models and map summaties as 
forecasts, rather than facts, appropriately cautious planning will likely ensue. 


One must also bear in mind data quality limitations that went in to the creation of the 
planning models. Digital terrain data is fairly good — 30 meter intervals between fairly 
accurate elevations — but vegetation data is rather poor. Vegetation data was derived in 
part from 500 meter grid cells of predicted natural vegetation. Thus, the worst common 
spatial denominator in the model is 500 meters. This has a major effect on the boundary 
between very different vegetation regimes, such as the piedmont to montane margin. 


The solution to many of these limitations lies in utilizing the model frequently. Actively 
noting inconsistencies (and consistencies) with forecast values will point out areas of 
poor baseline data, insufficient archaeological knowledge, or both. Both deficiencies can 
be remedied. Baseline data can be fixed on a local level, and more inventory in poorly- 
represented settings can be a management goal. From a land use perspective, confirming 
LOW forecast areas may be the highest priority. 


Maintaining the model is critical to its utility. Field protocols for gathering model data 
are straightforward. A simple tally sheet for each inventory can be created that 
summarizes the areal coverage in each model zone, and the revealed size density within 
each zone. Each inventory and resource should be held in GIS, verified, and flagged as 
not having contributed to the current generation of the model. Periodically, the model 
maintainers need to review new information and decide what effort should be put in to 
model revisions. This could be as simple as just changing the forecast maps without 
Statistical re-analysis or as comprehensive as running entirely new tallies and contrasts. 


Resource distributions, overall, are relatively sensible. It is not difficult to understand the 
distribution of historic resources within the sensitivity model. They tend to lie near to 
water and near to transportation routes. This generalization shows clearly within each of 
the analytical unit studies. Nevertheless, as a forecast of where significant or interesting 
historic resources will be found, the map layers should be used cautiously. For example, a 
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recent inventory of 10,000 acres within the Ely Field Office management area revealed 
few large historic sites, but dozens of dispersed small scatters of cans and some 
glassware. These were likely sheep camps. Careful analysis and locating of these 
seemingly insignificant, uninteresting, sites revealed a good deal about the settlement 
pattern of early twentieth century sheep-rearing. Each individual site would not have 
been considered significant; together they are a potential National Register landscape 
(Clay, persona) communication 2002). 


Prehistoric resources are more difficult to understand in simple factorial ways. The 
variation from one study unit to the next is somewhat unexpected. In some units, there 
are “sensible” reasons. For example, marshes and dunes (which lie along flat-piedmont 
interfaces) are important areas for food resources in the Bonneville Basin. Sites tend to be 
more frequent in these places. Other results are less “sensible”. Why should sagebrush 
flats be more likely to contain archaeology than piedmont in several of the Nevada study 
units? Why is water sometimes a negative factor? As promised, we offer no answers for 
these questions. They do make clear the importance of continuing to develop explanatory 
(causal, deductive) models alongside of correlation forecasts. The kind of study presented 
here, (an example of the latter activity) will be improved by creating forecasts from a 
better understanding of the rationale behind prehistoric behavior. 


Planning Models As Tools 


The utility and limitations for the planning model have been discussed above, but it is 
worth reiterating these again. The appropriate use of these tools is: 


Long range planning 

High probability relates to greatest likely overall expense 

Low probability equates with fewer resources, lower overall expense. 

If fewer sites are encountered, then testing, mitigation costs are reduced. 

Low probability does not mean no sites and does not obviate the need for 

fieldwork. But fieldwork should be faster and cheaper, on average. 

e As model is verified further, cultural resource managers may want to examine 
different level of investigation within low probability areas. 

© Models and forecasts articulate current state of knowledge. Thus they need 

maintenance. 


Closing Perspectives 


The GBRI model study was, we think, successful. Success is always a relative term. For 
this study, one goal was to evaluate the feasibility of building extensive landscape-level 
models. A feasibility evaluation was achieved: yes, this is feasible. Another goal was then 
to build such a model. This goal was also met in a series of analytical unit studies. A third 
goal was to examine how the model could best be used. This goal was partially met. 
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Above, we proposed ways in which the model and forecasts can be used and improved. 
This should be a continual work in progress as managers and scientists work together to 
improve ways in which decisions are made. 


Future work and elaboration needs to focus on improving survey methods, continually 
evaluating the model and forecast reliability, “fixing” the forecasts, and revising the 
model with better baseline environmental data. All of these future action 
recommendations pre-suppose agency use of the forecasts in the first place, and an 
audience for them. The agency staff, external researchers, and managers with whom we 
worked on this study are all sincerely interested in the success of landscape-level 
analyses. We think this is the best forecast, in itself, of continued use for this research. 
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